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PROPERTIES  OF  STEADY,  VISCOSITY-STRATIFIED 
FLOW  TO  A LINE  SINK 


CHAPTER  I 

INTRODUCTION 

Background 

A fluid  In  which  the  viscosity  varies  in  the  vertical  direction 
is  defined  as  being  "viscosity  stratified."  The  effects  of  viscosity 
stratification  are  investigated  in  this  report  using  a two-dimensional, 
non-diffuse,  laminar  flow  mathematical  model  which  was  developed  for  low 
Reynolds  numbers. 

Flows  in  which  density  stratification  exist  have  been  of  con- 
siderable engineering  interest.  Examples  of  these  flows  include  selec- 
tive withdrawal  from  reservoirs,  density  currents  in  reservoirs,  and 
salt  wedges  in  estuaries.  In  these  instances,  the  variation  in  vis- 
cosity of  the  fluid  has  been  ignored  in  mathematically  modeling  the 
flow  of  the  fluid  to  a sink  with  negligible  error. 

In  some  fluids  large  variations  of  the  rheological  properties 
of  the  fluid  with  depth  can  exist.  Walskl  and  Schroeder  (1978),  in 
attempting  to  model  the  behavior  of  dredged  material  flowing  over  a 
weir  at  a dredged  material  containment  area,  speculated  that  some  of 
the  differences  between  observed  and  predicted  velocity  profiles  were 
due  to  the  increased  viscosity  in  the  lower  layers  of  material  in  the 
containment  area.  The  ability  to  better  describe  the  behavior  of 
fluids  which  are  viscosity  stratified  could  result  in  more  efficient 
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design  of  weirs  at  these  containment  areas  or  in  any  separation  process 
in  which  viscosity  gradients  exist. 

The  problem  is  complicated  by  the  fact  that  viscosity  gradients 
are  usually  accompanied  by  density  gradients.  The  literature  does  not 
contain  solutions  of  the  equations  of  motion  for  two-dimensional  flow 
to  a sink  with  both  viscosity  and  density  gradients  existing.  There 
are  no  criteria  to  predict  when  viscosity  stratification  affects 
velocity  profiles,  nor  has  the  magnitude  of  these  effects  been 
quantified. 

Objectives 

The  objectives  of  this  study  are  to  experimentally  determine 
the  magnitude  of  the  effect  of  viscosity  stratification  on  two-dimen- 
sional flow  to  a line  sink  and  to  develop  a technique  for  mathemati- 
cally modeling  such  flow.  (Note  that  in  the  above  objective  a weir 
is  considered  a special  case  of  a line  sink,  i.e.,  surface  sink.)  This 
study  develops  not  only  criteria  for  predicting  when  viscosity  strati- 
fication is  significant  but  also  for  predicting  the  velocity  profiles 
in  those  cases. 

Overview  of  Report 

In  order  to  adequately  cover  this  topic,  a fairly  lengthy 
report  is  necessary.  The  chapters  have  been  written  so  that  they  may 
be  read  individually.  The  only  exceptions  to  this  are  that  Chapter  III 
(Basic  Equations  and  Concepts)  is  a prerequisite  for  Chapters  IV,  VI, 


VII,  VIII,  and  IX,  and  Chapter  V is  similarly  a prerequisite  for  Chapters 
VI  through  IX. 

The  literature  review  outlines  state-of-the-art  at  the  time  this 
study  was  done.  Chapter  III,  Basic  Equations  and  Concepts,  begins  with 
a list  of  assumptions  used  in  the  mathematical  model  development.  The 
horizontal  and  vertical  momentum  equations,  continuity  equation,  and 
species  continuity  equation  are  then  simplified  and  nondimensionalized 
to  yield  the  equations  to  be  solved  in  later  parts  of  the  report.  The 
boundary  conditions  and  numerical  techniques  used  are  also  discussed. 

One-dimensional  models  are  developed  in  Chapter  IV  for  a variety 
of  viscosity  gradients  and  boundary  conditions  to  illustrate  the  effect 
of  viscosity  stratification.  Skewness  of  velocity  profiles  is  defined 
as  an  indicator  of  the  effect  of  viscosity  stratification  on  velocity 
profiles  in  Chapter  V. 

Chapter  VI,  Linear  Viscosity  Gradients  Neglecting  Buoyancy 
Forces,  contains  solutions  to  the  equations  developed  in  Chapter  III 
first  considering  only  viscous  forces,  then,  later  inertial  forces. 

The  equations  are  first  solved  for  duct  boundary  conditions,  then  for 
boundary  conditions  corresponding  to  a fixed  upstream  wall. 

In  Chapter  VII,  Model  for  a Linear  Gradient  with  Buoyancy 
Forces,  the  equations  are  solved  with  the  buoyancy  force  term  included. 
Criteria  for  the  existence  of  viscosity  stratification  are  developed. 
Chapter  VII I addresses  the  same  problem  as  the  previous  chapter  except 
that  the  viscosity  gradient  Increases  logarithmically  with  depth.  The 
concept  of  "incipient  skewness"  is  Introduced  to  describe  cases  in 
which  buoyancy  forces  predominate  near  the  top  of  the  fluid  while 
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viscous  forces  predominate  below. 

Chapter  IX  contains  a discussion  of  the  effect  of  some  of  the 
simplifications  made  on  the  previous  solutions.  These  include  the 
effects  of  the  uniform  upstream  velocity  assumption,  the  fixed  back 
wall,  and  neglecting  inertial  forces. 

Laboratory  procedures  are  presented  in  Chapter  X,  including  a 
discussion  of  equipment,  fluids,  and  procedures.  Experiment  results 
are  then  presented  in  Chapter  XI.  First,  data  for  the  runs  are  given. 
This  is  followed  by  a discussion  of  the  agreement  between  the  mathe- 
matical model  and  experimental  results  and  an  analysis  of  errors. 

In  Chapter  XII,  the  model  is  then  applied  to  weirs.  The  effect 
of  viscosity  stratification  on  concentrations  of  a fluid  flowing  over  a 
weir  Is  discussed. 

The  properties  of  non-Newtonian  flow,  in  particular  a Bingham 
plastic  with  yield  stress  increasing  with  depth,  are  presented  in 
Chapter  XIII.  The  implications  of  this  on  design  are  discussed  for 
yield  stress  increasing  linearly  and  logarithmically  with  depth. 
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CHAPTER  II 


LITERATURE  REVIEW 

Existing  literature  in  several  areas  of  fluid  mechanics  was 
surveyed  to  locate  solutions  to  the  problems  of  selective  withdrawal 
from  a line  sink  in  a two-dimensional,  continuously  stratified  fluid 
with  varying  viscosity.  While  no  papers  could  be  found  which  directly 
addressed  the  problem,  a large  number  were  available  on  areas  closely 
related  to  the  problem.  These  papers  are  grouped  into  the  eight 
categories  presented  below. 

General  Reviews 

One  of  the  earliest  reviews  on  stratified  flow  through  the 
late  1950's  was  prepared  by  Harleman  (1961).  The  selective  withdrawal 
work  at  this  time  was  limited  generally  to  two-layer  systems. 

Yih  (1965,  1967,  1970)  published  a number  of  works  on  strati- 
fied flow.  In  general,  Yih's  work  on  selective  withdrawal  is  limited 
to  the  inviscld  flows  with  a linear  density  gradient  far  upstream. 

Brooks  and  Koh  (1968)  presented  an  excellent  overview  of  the 
state-of-the-art  of  selective  withdrawal  models  in  the  late  1960's. 
They  dealt  chiefly  with  reservoir  stratification,  and  did  not  consider 
variations  in  viscosity. 

Turner's  (1973)  text.  Buoyancy  Effects  in  Fluids,  does  not 
dwell  on  the  subject  of  selective  withdrawal,  but  it  does  Include  a 
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discussion  of  the  differences  between  solutions  ignoring  diffusion, 
which  rely  on  the  Richardson  and  Reynolds  number,  and  those  that 
account  for  diffusion,  which  rely  on  a Rayleigh  number. 

Early  Research 

Most  early  research  on  stratified  fluids  was  concerned  with 
the  problem  of  stability,  e.g.,  Taylor  (1931).  When  viscosity  was  con- 
sidered as  in  Keulegan  (1949),  it  was  treated  as  being  constant. 

The  problem  of  selective  withdrawal  was  first  approached  by 
Gariel  (1949)  and  Craya  (1949).  Craya  pointed  out  that  viscous  forces 
will  predominate  as  one  moves  away  from  the  sink.  Gariel  was  first  to 
show  the  shape  of  the  velocity  profile  in  this  zone. 

Harleman,  Gooch,  and  Ippen  (1958)  studied  withdrawal  from  a 
sluice  gate  and  found  the  withdrawal  layer  thickness  to  be  dependent 
on  the  densimetric  Froude  number.  Harleman,  Morgan,  and  Purple  (1959) 
then  extended  this  work  to  different  types  of  withdrawal  structures. 

Huber  (1960)  produced  a solution  for  a two-layer  fluid  flowing 
through  a line  sink  using  a relaxation  technique  and  found  the  flow  to 
be  governed  by  the  densimetric  Froude  number. 

Martin  and  Schuster  (1956)  produced  the  first  study  of  strati- 
fied flow  over  a weir.  They  were  the  first  to  speculate:  "The 
effect  of  increasing  viscosity  with  depth,  although  not  a predominant 
factor  of  withdrawal,  would  be  to  lessen  the  contribution  of  the 
lower  strata  to  the  total  outflow."  They  modeled  thermal  stratifica- 
tion with  salt  solutions.  They  speculated  on  using  glycerol  to  model 
the  viscosity  but  abandoned  the  idea  since  they  could  not  achieve  the 
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desired  density  and  viscosity  gradients  simultaneously.  They  therefore 
ignored  viscosity  in  the  analysis. 

Schlag  (1965)  did  the  most  thorough  study  of  two-layer  flow  over 
a weir.  He  was  first  to  study  the  case  in  which  the  weir  did  not  extend 
completely  across  the  channel.  He  produced  a set  of  nomograms  to  de- 
termine withdrawal  depth. 

The  above  studies  were  concerned  primarily  with  predicting  the 
thickness  of  the  withdrawal  layer.  No  attempt  was  made  to  predict 
velocity  profiles  except  to  indicate  average  velocity  in  each  layer  of 
a two-layer  fluid. 


Empirical  Studies 

In  a series  of  studies  produced  at  the  U.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES),  Grace  (1971)  and  Bohan  and  Grace  (1973) 
produced  an  empirical  model  which  predicted  both  depth  of  withdrawal 
layer  and  velocity  profiles  for  flow  over  sharp-  and  broad-crested  weirs. 
The  model  has  been  tested  for  a wide  range  of  cases  both  in  the  laboratory 
and,  to  a lesser  extent,  in  the  field. 

The  model  was  later  used  to  model  turbidity  in  an  impoundment 
by  Fontane,  Bohan,  and  Grace  (1973).  This  application  did  not  consider 
the  effect  of  turbidity  on  viscosity,  and  the  model  does  not  account 
for  the  effect  of  viscosity  on  the  velocity  profile. 

Johnson  (1974)  produced  a model  very  similar  to  that  of  Bohan 
and  Grace.  It  also  neglected  viscous  forces. 

Walski  and  Schroeder  (1978)  applied  the  WES  selective  withdrawal 
to  flow  over  a weir  in  a dredged  material  containment  area.  They  found 


the  model  after  some  empirical  corrections  could  describe  the  observed 
velocity  profiles  measured  at  several  containment  areas.  They  specu- 
lated that  the  corrections  to  the  model  were  due  to  the  fact  that  the 
suspended  solids  did  not  only  create  a density  gradient  but  a signif- 
icant viscosity  gradient. 

Several  models  of  flow  in  sedimentation  tanks  have  been  pre- 
pared by  Ingersoll,  McKee,  and  Brooks  (1956)  and  Clements  (1966).  These 
did  not  consider  the  effects  of  stratified  flow. 

Inviscid  Flow  Models 

Debler  (1959),  working  with  Yih  (1958),  showed  that  selective 
withdrawal  would  occur  for  densimetric  Froude  numbers  less  than  1/ir. 
Debler' s technique  for  filling  the  flume  served  as  a model  for  later 
researchers,  including  this  study. 

Yih's  solution  was  invalid  for  values  of  the  densimetric  Froude 
number  less  than  1/ir.  Kao  (1965)  was  able  to  solve  the  problem  using 
a free  streamline  solution  to  model  this  case.  His  theoretical  results 
roughly  agreed  with  Debler 's.  The  differences  could  be  explained  by 
the  fact  that  Kao  Ignored  viscous  effects. 

Trustrum  (1964)  showed  the  relationship  between  blocked  flows 
in  stratified  and  rotating  fluids.  Later  Janowitz  (1973)  showed  that 
introduction  of  viscosity  could  have  a significant  impact  on  flow  over 
a vertical  barrier  in  an  unbounded  case. 

In  a series  of  studies  at  the  University  of  Canterbury,  Wood 
(1968),  Wood  and  Lai  (1972a,  1972b),  and  Lust  and  Wood  (1974)  solved  a 
two-  (or  multi-)  layer  energy  equation  for  flow  through  a gradually 
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tapered  weir,  sluice  gate  and  orifice.  They  determined  the  discharge 
from  each  layer. 

Armi  (1975)  solved  the  inviscid  energy  and  momentum  equations 
for  flows  through  contractions  and  over  broad-crested  weirs  for  the  two- 
layer  case.  The  solutions  were  dependent  on  the  Froude  numbers.  He 
ignored  viscous  forces. 

Viscous  Flow  Models 

Koh  (1964,  1966a)  was  first  to  solve  the  problem  of  viscous, 
diffuse  flow  toward  a line  sink.  He  used  Boussinesq  and  boundary  layer 
assumptions  and  the  slow  flow  assumption  to  simplify  his  equation. 
(Boussinesq  assumption  states  pQ  + Ap  ^ pq*  if  stratification  is  small.) 
Similar  assumptions  are  used  in  this  study.  Koh  then  was  able  to  define 
a similarity  transformation  to  reduce  the  equation  which  he  solved  using 
a perturbation  scheme. 

His  solution  showed  that  the  flow  would  be  symmetric  about  the 
elevation  of  the  sink.  His  solution  predicted  a large  forward  flowing 
zone  with  a small  backflow. 

Koh  made  the  assumption  that  the  stratification  is  small  which  per- 
mitted neglecting  derivatives  of  density  with  depth  in  inertial  terms.  At 
the  same  point  in  his  development  he  eliminated  derivatives  of  viscosity 
with  depth,  using  the  implicit  argument  that  since  the  density  stratifica- 
tion is  small,  the  viscosity  stratification  must  also  be  small.  While 
this  assumption  is  valid  for  some  fluids,  especially  salt  stratification 
which  Koh  and  others  used,  it  does  not  follow  that  it  will  be  true  for  all 
fluids. 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  E). 


10 


The  same  reasoning  has  been  used  by  all  later  researchers. 
Among  the  advances  presented  in  this  study  are  criteria  stating  when 
this  assumption  is  no  longer  true. 

In  contrast  to  most  other  researchers,  Koh  carried  the  term 


where 

p * viscosity 
y = vertical  distance 
u = horizontal  velocity 

through  most  of  his  development  before  eliminating  it  with  what  is 

essentially  a Boussinesq  assumption  for  viscosity.  (Most  other  re- 

2 

searchers  begin  their  development  using  pV  u to  represent  the  viscous 
stress  tensor.)  Koh  writes  this  assumption  as 

| Ap | <<P0  (2) 

where 

Ap  ■ change  in  viscosity  due  to  stratification 
Pq  ■ characteristic  viscosity 

For  a fluid  with  a linear  temperature  gradient  from  4 to  20°C, 

PQ  is  1.235  cp  and  Ap  is  0.565  cp.  Koh's  assumption,  therefore,  is 
not  always  true.  The  results  are  even  more  dramatic  when  stratified 
fluids  with  widely  varying  viscosities  are  used. 

The  techniques  for  measuring  velocity  developed  by  Koh  for  this 


study  have  been  used  by  a large  number  of  later  researchers,  including 
this  study,  with  some  modification.  It  consisted  of  dropping  dye 
pellets  into  the  fluid  and  photographing  them  at  later  times.  The 
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displacement  of  the  streaks  indicates  the  velocity. 

Gelhar  and  Mascolo's  (1966)  work  in  selective  withdrawal,  and 
Long's  (1962)  work  in  "jets,"  used  assumptions  similar  to  Koh  (1964) 
except  that  they  ignored  the  effect  of  diffusion  which  in  most  cases 
is  small.  They  found  velocity  profiles  which  were  virtually  identical 
to  Koh's,  which  indicates  that  the  effect  of  diffusion  was  probably 
small  in  Koh's  experiment. 

Gelhar  and  Mascolo  developed  the  differential  equation: 


i-X  + I IX  o 

3Y4  V 3X  ' 

where 


(3) 


S'  = stream  function 

X,Y  = horizontal  and  vertical  distances,  respectively 

f = Z£iP 
p 3 S' 
o 

p = density 

v = kinematic  viscosity 

which,  except  for  the  fact  that  it  assumes  constant  viscosity,  is  the 
equation  to  be  solved  in  this  study. 

Imberger  (1970)  and  Imberger  and  Fischer  (1970)  showed  that 
several  zones  governed  by  different  force  balances  exist  in  stratified 
flow  to  a sink.  While  they  defined  several  zones,  the  two  most  criti- 
cal are  the  outer  zone  in  which  a viscous-buoyancy  balance  predominates, 
and  an  inner  zone  in  which  the  inviscld  inertia-buoyancy  balance  pre- 
dominates. 

Unlike  Koh  (1964,  1966a)  and  Gelhar  and  Mascolo  (1966), 
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where  D = molecular  diffusion  coefficient. 


Imberger  (1970)  used  a finite  duct  model  instead  of  a semi- inf inite 
fluid  model,  although  he  did  allow  slip  along  the  duct  wall.  He  found 
the  most  important  parameters  to  be  the  Raleigh  number  ^VvD/egi2  ) and 
the  Schmidt  number  (Vv/D  ),  where  D = molecular  diffusion  coefficient. 

While  he  was  able  to  solve  the  outer  equations  using  a Fourier 
series  technique,  he  needed  an  approximate  integral  technique  for  the 
inner  layer.  His  laboratory  results  showed  excellent  agreement  with 
the  integral  solution. 

Imberger  and  Fischer  (1970)  criticize  Gelhar  and  Mascolo  (1966) 
for  ignoring  diffusion,  especially  since  they  were  solving  the  problem 
in  an  infinite  domain.  They  criticize  Koh's  (1964,  1966a)  solution 
near  the  sink  since  it  ignores  inertia  and  far  from  the  sink  since 
the  linear  density  profile  is  inconsistent  with  an  infinite  domain. 

In  a series  of  papers  prepared  at  the  University  of  Wisconsin, 
Walesh  and  Monkmeyer  (1973),  He  (1973),  and  Clark  (1974)  studied  slow, 
viscous,  stratified  flow  to  line  sinks  for  several  reservoir  geometries. 
Their  solution  method  was  similar  to  Gelhar  and  Mascolo  (1966)  although 
they  did  extend  the  method  to  axisymmetric  cases  in  cylindrical  co- 
ordinates. They  also  analyzed  shear  stress  in  stratified  flow.  They 
developed  an  empirical  relationship  for  the  f term  in  Gelhar  and  Mas- 
colo's  model.  Their  model  has  a realistic  boundary  condition  of  no 
flow  at  solid  boundaries  although  they  treat  viscosity  as  a constant. 

King,  Norton,  and  Orlob  (1973)  developed  a finite  element  model 


which  solves  the  Reynolds  equation  for  variable  stratification  over  a 
broad-crested  weir.  They  replaced  viscosity  with  eddy  viscosity  which 
they  determined  empirically.  They  only  reported  the  results  for  two 
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cases,  homogeneous  flow  and  a two-layer  system  with  a large  density 
gradient.  The  model  agreed  reasonably  well  with  the  data  provided  by 
Bohan  and  Grace.  Since  the  solution  was  developed  for  high  Reynolds 
numbers,  it  was  more  closely  related  to  the  inviscid  case  than  the  vis- 
cous. The  authors  reported  difficulty  in  applying  the  model  for  low 
Reynolds  numbers. 

One  source  which  contains  some  examples  of  steady  one-dimen- 
sional flow  in  the  form  of  example  problems  is  Bird,  Stewart,  and  Light- 
foot's  (1960)  Transport  Phenomena . While  it  does  not  contain  any  two- 
dimensional  problems  with  varying  viscosity,  its  example  served  as  a 
guide  to  one-dimensional  solutions  presented  later  in  this  report. 

For  the  case  of  steady  seepage,  with  Peclet  number  (Jlu/D) 
greater  than  one,  Yih  (1961)  showed  that  the  problem  of  varying  vis- 
cosity can  be  solved  in  a straightforward  manner.  The  technique  in- 
volved solving  the  problem  of  potential  flow  in  a porous  media  for  a 
homogeneous  fluid.  The  velocity  determined  by  this  solution,  which 
Yih  called  the  velocity  of  associated  flow,  can  then  be  transformed  to 
the  velocity  at  any  point  by  multiplying  by  the  ratio  of  the  reference 
viscosity  to  the  viscosity  at  that  point.  This  approach  could  also 
apply  to  flows  in  Hele-Shaw  cells.  In  this  study,  the  derivatives 
in  the  direction  of  flow  and  vertical  direction  are  at  least  the  size 
of  derivatives  normal  to  the  flow  so  the  above  transformation  is  of 
little  use. 

Unsteady  Flow  Models 

Koh  (1966b),  and  later  Imberger  and  Fischer  (1970),  solved  the 
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initial  value  problem  for  the  inviscid  case.  By  matching  the  solution 
to  a steady-state  solution,  they  determined  the  time  for  a withdrawal 
layer  to  form.  This  result  is  used  in  this  study  to  check  to  see  if 
steady-state  conditions  were  reached  in  the  model  tests. 

In  a series  of  articles,  Kao,  Pao,  and  Wei  (1972),  Pao  and  Kao 
(1974),  Kao,  Pao,  and  Wei  (1974),  and  Kao  (1976)  developed  a numerical 
model  to  more  closely  study  the  mechanisms  involved  in  the  establish- 
ment of  stratified  flow  in  an  infinite  duct.  They  showed  that  steady- 
state  develops  as  a series  of  columnar  disturbances  propagating  out 
from  the  sink.  They  found  viscous  forces  to  be  unimportant  for  Reynolds 
numbers  greater  than  1000. 

Fandry  (1974),  Imberger  and  Fandry  (1975),  and  Imberger, 

Thompson,  and  Fandry  (1976)  extended  unsteady-state  solutions  from  the 
infinite  duct  problem  to  a tank  with  a back  wall.  They  showed  that 
reflection  of  the  disturbances  from  the  back  wall  can  result  in  a 
secondary  collapse  of  the  withdrawal  layer  depending  on  the  Prandtl 
number  (v/D)  and  a parameter  relating  inertial  and  buoyancy  forces 
(q(eg*2v2)-1^3) . 

Numerical  Methods 

This  study  does  not  involve  any  major,  new  developments  in 
numerical  solution  of  fluid  mechanics  problems.  The  basic  reference 
used  was  Carnahan,  Luther,  and  Wilkes  (1969),  Applied  Numerical  Methods. 
Another  reference  used  was  a state-of-the-art  paper  on  simulation  of 
viscous  flows  by  Orszag  and  Israeli  (1974) . 

Imberger,  Thompson,  and  Fandry  (1976)  point  out  that  for  flows  in 
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which  large  density  differences  develop,  an  Eulerian  solution  is  not  ap- 
propriate. They  used  a mixed  Eulerian-Lagrangian  solution  in  their  study 
similar  to  "cloud  in  cell"  techniques  in  which  particles  of  constant  den- 
sity are  tracked  through  an  Eulerian  field.  They  implied  that  this  tech- 
nique was  expensive. 


Rheological  References 

Significant  viscosity  gradients  in  fluids  can  be  achieved  by 
the  presence  of  colloidal  suspensions.  There  is  a linear  relationship 
between  concentration  and  viscosity  for  dilute  suspensions  as  given  by 
Einstein's  Law  of  Viscosity.  This  relationship  becomes  non-linear  at 
higher  concentrations  and  non-Newtonian  effects  predominate  at  even 
higher  concentrations.  Reviews  of  these  models  of  the  effect  of  sus- 
pensions on  viscosity  have  been  prepared  by  Frisch  and  Simha  (1956)  and 
Jeffrey  and  Acrivos  (1976). 

Sediments  were  found  to  have  Bingham  plastic  behavior  at  con- 
centrations greater  than  30  g /l  by  Krone  (1963a,  1963b)  and  Migniot 
(1968).  Krone  found  the  yield  stress  to  be  a function  of  concentration 
to  the  2.5  power  while  Migniot  determined  it  to  be  to  the  fourth  power. 

The  basic  references  used  for  determining  the  rheological  prop- 
erties of  model  fluids  used  in  this  study  were  Metzner  (1956),  Reiner 
(1960),  and  Van  Wazer  et  al.  (1963). 

Summary 

No  solution  could  be  found  to  the  problem  of  tvo-dlmensional, 
stratified,  non-Newtonian  flow  to  a line  sink. 
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In  general,  the  effects  of  a viscosity  gradient  on  the  velocity 
distribution  in  stratified  flow  have  been  ignored,  since  in  the  cases 
studied  thus  far  the  effects  of  a density  gradient  are  more  dramatic 
than  those  caused  by  the  viscosity  gradient.  In  particular,  in  the 
case  of  slow,  viscous  flow  to  a line  sink,  previous  researchers  have 
made  the  assumption  that  the  density  stratification  is  small,  which 
allows  them  to  invoke  the  Bousslnesq  assumption  (pq  *v  PQ  + Ap). 

They  have  then  concluded  that  since  the  density  stratification 
is  small,  then  the  viscosity  stratification  must  also  be  small 
(y  *v  yQ  + Ay).  In  this  report,  criteria  are  developed  to  determine 
when  the  assumption  of  small  viscosity  stratification  is  not  true  and 
techniques  for  predicting  the  velocity  distribution  are  presented  for 


that  case. 


I 

CHAPTER  III 

BASIC  EQUATIONS  AND  CONCEPTS 

In  this  chapter,  the  scope  of  the  study  is  precisely  defined  by 
listing  the  assumptions  used  in  the  study.  The  equation  of  motion  to  be 
solved  is  presented  with  special  emphasis  on  the  Dt^/3 y term  in  the 
horizontal  momentum  equation  which  accounts  for  viscosity  stratifica- 
tion (t  is  viscous  stress  tensor).  Notational  conventions  and  dimen- 
sionless groups  are  discussed  and  the  equations  are  presented  in  a 
manner  so  that  it  is  unnecessary  to  repeat  these  preliminaries  in 
later  chapters  of  the  report. 

Assumption 

The  following  assumptions  are  used  to  define  the  scope  of  the 
study.  The  first  set  of  assumptions  are  true  in  a large  number  of 
practical  cases  for  laminar  flow. 

1.  Flow  is  laminar  (no  turbulent  momentum  transfer  terms  in 
equation) ; 

2.  Stratification  is  Boussinesq  (pq  'v  pq  + Ap); 

3.  Flow  is  non-diffuse  (D  ^ 0); 

4.  Flow  is  steady  in  duct  models  and  quasi-steady  in  tank  models 
(derivatives  with  respect  to  time  are  neglected) ; 

5.  Flow  is  incompressible  (V  • il  ■ 0) ; 

6.  Inertial  terms  in  vertical  motion  equation  are  neglected  (ac- 
celerations are  only  important  in  horizontal  direction) ; 
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7.  Density  is  a linear  function  of  concentration  of  species 
affecting  density  (density  can  be  used  in  place  of  concentra- 
tion in  species  continuity  equation); 

8.  There  are  no  body  forces  in  horizontal  direction  (g^  = 0); 

9.  Viscosity  is  independent  of  pressure,  (8p/3p  = 0); 

10.  Flow  is  two-dimensional  (derivatives  in  the  lateral  direction 
can  be  ignored) . 

A second  set  of  assumptions  may  not  generally  be  true  in 
practical  problems.  Nevertheless,  the  assumptions  are  useful  in  making 
the  mathematics  tractable.  Each  of  these  assumptions  will  be  relaxed 
during  the  study  to  determine  its  effects. 

1.  Flow  is  Newtonian  (Newton's  Law  of.  Viscosity  can  be  applied); 

2.  Flow  is  slow  (u  • Vu  = 0) ; 

3.  Density  gradient  is  linear  (3p  / 3y  = sp  = constant); 

s o 

4.  Line  sink  is  located  at  mid-depth  in  the  model  (horizontal 
boundaries  have  equal  impact  on  solution) ; 

5.  Boundary  layer  type  assumptions  can  be  used  (3u/3y  » 3u/3x). 

Development  of  Basic  Equations 

The  following  equations  are  a direct  result  of  the  first  set  of 
assumptions.  The  terms  are  defined  in  Figure  1. 
x-momentum 


(5) 
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continuity 


3u  3v 
3x  3y 


(6) 


where 


species  (density)  continuity 

ui£  + v ip+!!£  - c 

3x  3y  3y 


pT  ■ p0  + p8(y)  + p<x»y> 

po  **  density  upstream  at  sink  elevation 


(7) 


3p 


constant 


PT  - Pg(y)  + p(x,y) 
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The  above  notational  conventions  have  the  effect  of  considering 
density  and  pressure  only  as  they  are  affected  by  motion,  i.e.,  sul 
tracting  out  static  pressure  gradient.  The  above  equations  are  essen- 
tially the  same  as  those  used  by  Imberger  (1970)  with  the  exception  that 
viscosity  is  now  allowed  to  vary  and  time  derivatives  are  neglected  in 
this  study. 

In  the  past  studies  the  term  has  been  treated  as  the 
term  and  replaced  by 


yx 


xx 


3t 

...g 

3y 


-v 


32u 

3y2 


(8) 


based  on  Newton's  Law  of  Viscosity  for  fluids  with  constant  viscosity. 
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To  allow  for  varying  viscosity 
3t 


3y 


3_ 

8y 


-y 


fau  . 3v  |1 

[ay  + 3xJ] 


2 

3 u 

ay 


2 

_3_y  3u  jly  3v  _ 3 v 

3y  3y  3x  3x  ^ 3x3y 


(9) 


The  lat>c  two  terms  can  be  eliminated  with  a boundary  layer  type  assump- 
tion. It  is  because  of  the  term  3y/3y  (3u/3y)  and  the  fact  that  y is 
2 2 

not  constant  in  y 3 u/3y  that  viscosity  stratification  can  affect  flow. 

This  study  will  examine  the  effect  of  retaining  the  3y/3y  term, 
which  will  be  called  the  viscosity  stratification  term.  Two  special 
cases  will  be  studied  extensively.  The  first  is  a linear  viscosity 
gradient 

V “ y (10) 


where  yQ  = viscosity  at  sink  elevation. 

The  second  is  a logarithmic  increase  of  viscosity  with  depth 
y = y e 1J 


(ID 


- slope  of  log  viscosity  curve. 


where  c^ 

Allowing  for  this  viscosity  stratification,  the  x-momentum 
equation  becomes 


(• 


3u 

3x 


+ v 


+ !* 

dx 


(32u  . 32u  1 . 3y  3u 

3x2  3y2J  * * 


(12) 


2 2 

Note  that  the  3 u/3x  term  will  be  eliminated  from  most  equations  by 
the  boundary  layer  assumption. 


Non-Dimenslonallzlng 


Throughout  this  report  the  equations  are  solved  in  their 


A 


< 
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dimensionless  form.  The  equations  are  non-dimenslonalized  using  a 

characteristic  two-dimensional  flow  at  the  sink  (q) , length  (1),  and 

2 

viscosity  (pq)  with  units  L /T,  L,  and  M/LT,  respectively.  When  the 
substitutions 

U = ui/q,  V - vi/q,  X *=  j,  Y - ^ 


are  made,  equations  4,  6 and  7 become: 
x-momentum 

qPo  fIT  3U  . „ 3U  ] . i3  9p  y / 32U  . 32uV  l 3p  3U 
Po  ^ 3X  3Y  J qyQ  3x  ' yo  ^2  3y  3Y 

continuity 


(13) 


3U  3V 
3Y  3Y 


0 


(14) 


species  continuity 


The  y-momentum  equation  is  unchanged. 

Note  that  capital  letters  will  be  used  to  designate  the  dimen- 
sionless value  of  the  variables.  No  attempt  was  made  above  to  non- 
dlmensionalize  the  pressure  and  buoyancy  terms  since  these  will  be 
handled  differently  in  the  various  cases. 

The  above  equations  will  be  used  as  a starting  point  for  many 
of  the  following  chapters.  Note  that  the  equations  contain  the  Reynolds 

number  (qp  /y  ).  In  most  cases  this  will  be  considered  small  and  can 
o o 

be  eliminated  from  the  equations. 

The  most  important  parameter  for  the  purposes  of  this  study  is 
l_  3u 

y0  »y 


which  is  an  indicator  of  the  strength  of  the  viscosity  stratification. 
It  will  give  rise  to  two  dimensionless  groups  depending  on  whether 
3u/3y  is  constant,  i.e.,  gradient  is  linear. 
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Simplification  of  Equations 


To  solve  the  system  of  equations,  a single  equation  with  one 
unknown  (V)  was  constructed  from  the  four  coupled  equations  (5,  13,  14, 
15)  with  four  unknowns  (U,  V,  p,  P) . This  is  done  below  for  the  case 
in  which  the  Reynolds  number  is  small  and 


2 2 
O « jTU 

2 2 
3X  3Y 


(16) 


First,  pressure  can  be  eliminated  from  the  momentum  equations  by 
differentiating  the  x— momentum  equation  by  y and  the  y— momentum  equation 
by  x to  give 

i3  32p  _ p 33U  + 2 3yi  32U  1 32y  3U  <t3g  3p  (1?) 

qpo  3x3y  3Y3  % 3y  ^2  y<>  3y2  3Y  " ‘ qvo  3x 

The  continuity  equation  implies  the  existence  of  a stream 
function  y such  that 


sy  3y 

3y  ” ’ 3x 


-V 


(18) 


Substituting  this  into  the  above  equation  and  rearranging  terms  gives 

+ (19) 

3Y4  v 3y  3Y3  v 3y2  3Y2  qv  3x 

The  species  continuity  equation  can  be  rearranged  (for  U f 0) 


as 


U “ ■ u h (p  + ps} 


(20) 
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While  3/3y  (p  + p ) is  not  known,  it  is  known  that  the  fluid  is 


non-diffuse  (density  constant  on  a streamline)  so  that 


3 (P  + Pg) 
3Y 


O T O / , N 

37  “ 37  (p  + Ps1* 


(21) 


At  some  point  very  far  from  the  sink,  the  density  field  is  not  affected 
by  flow  and  it  is  possible  to  write 


3p  3 (p  + p ) 

3 S 

3Y  " 


3Y 


11  = A 111 

3y  ~ i l 3Y 


(22) 


3 (P  + Ps) 


3Y 


is  a constant,  and 


3 (p 


1^1-  if”) 

1 l 3y  I 


« I 87  (p  + P8> 


Jlep 


o 3Y 


(23) 


where 


. . 1_  dp 
po  dy 

Equation  23  can  be  rearranged  (noting  3f/3Y  “ U)  and  inserted  into 
equation  21  to  give 

3 "epoU 

*7  (p  + p8)  ‘ ir 

J 00 

Substituting  back  into  equation  20  gives 

3p  . V 

3x  “ U epo 
00 

Substituting  for  3p/3x  in  the  motion  equation  gives 


(24) 


(25) 


h*2  A »iL  I??  ♦ ii  '»  »’ 

3YA 


t Ep^g  U, 

SL  - 0 


V 3Y3  V 3y2  3Y2  qPoU-  W *X 


(26) 


The  above  equation  is  the  model  for  stratified  flow  that  is  used 
in  the  remainder  of  the  report  whenever  buoyancy  forces  are  important. 


This  equation  is  used  as  the  starting  point  for  the  models  developed 
in  Chapters  VII,  VIII,  IX,  XI,  and  XII. 

The  solutions  are  shown  not  to  be  very  sensitive  to  the  term 
pU(o/yo  since  it  deals  with  upstream  boundary  conditions  which  are  of 
only  minor  importance  since  the  equation  is  solved  as  an  initial  value 
problem  in  x.  U is  usually  set  to  one  to  simplify  the  solution  pro- 
cedure . 

Similarly,  y/yQ  is  replaced  by  the  appropriate  y(Y)  for  that 
model.  Note  that  since  the  flow  is  non-diffuse,  y is  constant  on  a 
streamline  and  is  a function  of  V not  Y.  Using  y as  a function  of  Y 
implies  that  the  viscosity  field  is  fixed.  The  effect  of  this  assump- 
tion is  small  and  is  discussed  in  Chapter  IX. 

Boundary  Conditions 

Each  of  the  cases  studied  in  this  report  requires  a slightly 
different  set  of  boundary  or  initial  conditions.  These  conditions  are 
categorized  in  this  chapter  and  are  developed  in  greater  detail  as  each 
case  is  analyzed. 

In  the  one-dimensional  models  the  only  question  is  whether  the 
fluid  is  allowed  to  slip  at  the  boundaries.  In  general,  it  will  not. 

In  two-dimensional  flow  problems,  the  distinction  between  duct  and  tank 
models  must  be  made.  The  flow  originates  at  x » « for  the  duct 
model , while  for  tank  models  there  is  a fixed  back  wall  and  the  flow 
originates  at  y » 1/2  (Y  = 1/2).  These  ca.es  are  shown  in  Figure  2. 

A much  more  important  question,  especially  in  the  duct  model, 
is  whether  the  fluid  is  allowed  to  slip  along  the  wall.  Most  earlier 
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Fig.  2.  Boundary  conditions 


w- 
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researchers  were  interested  in  withdrawal  layers  much  smaller  than  the 
height  of  the  duct  and  hence  allowed  the  flow  to  slip.  For  viscosity 
stratification  to  become  important,  it  is  important  that  the  fluid  not 
be  allowed  to  slip. 

In  general,  the  problem  is  a boundary  value  problem  in  the 
y-direction.  In  the  x-direction,  it  is  a boundary  value  problem  in 
the  tank  model  and  an  initial  value  problem  in  the  duct  model.  The 
boundary  conditions  for  the  duct  model  are  given  in  Figure  3. 


Numerical  Techniques 


Numerical  methods  for  solving  differential  equations  are  re- 
quired for  a number  of  problems  in  this  study.  Finite  difference 
methods  are  used  to  reduce  the  problem  to  a system  of  algebraic  equa- 
tions which  are  solved  by  the  Gauss-Seidel  Iterative  Method. 

The  numbering  convention  for  the  grids  is  shown  in  Figure  4. 
The  finite  difference  approximations  for  derivatives  used  in  this 
study  are  listed  below  using  U as  the  dependent  variable  and  Y as  the 
independent  variable. 

First  derivative  (boundary  value  problem) 


au  Vi-Vi 

3Y  24  Y 

First  derivative  (initial  value  problem) 
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Second  derivative 


- 2U1  * Vi 

(AY)2 


Third  derivative 


Vl  - 3U1  * 3U1-1  - U1-2 

(AY)3 


Fourth  derivative 


(29) 


(30) 


V2  - *Vl  + 6U1  - 4U1-1  + U1-2 

(AY)4 


(31) 


In  problems  that  were  solved  in  terms  of  U,  no  slip  boundary 

conditions  were  handled  by  setting  U = 0 at  the  boundary.  In  problems 

solved  in  terms  of  Y,  the  stream  function,  this  condition  became 

3Y/3Y  » 0.  This  implied  YN_^  could  then  be  substituted  for  Y^t  in 

4 4 

formulas  for  higher  derivatives  (e.g.,  3 Y/3Y  ) for  a boundary  at  j = N. 

Slip  boundary  conditions  in  a Newtonian  fluid  can  be  interpreted 
2 2 

as  3U/3Y  = 0 or  3 Y/3Y  “0.  A substitution  as  given  in  the  preceeding 

paragraph  suffices  for  the  condition  in  terms  of  U.  For  the  case  of  Y, 

4 4 

the  higher  derivative  (e.g.,  3 Y/3Y  ) is  written 


and  32Y/3Y2  is  substituted  in  the  32/3Y2  formula  for  j - N. 

All  numerical  solutions  were  programmed  in  FORTRAN  IV  and 
solved  on  a UNIVAC  1108  computer.  Descriptions  and  flow  chart  of  the 
programs  are  given  in  Appendix  A. 


CHAPTER  IV 


ONE-DIMENSIONAL  MODELS 


Velocity  profiles  are  developed  in  this  part  of  the  report  for 
some  simple  one-dimensional  flows.  The  profiles  developed  here  are 
used  to  demonstrate  the  effect  of  viscosity  gradients  on  velocity,  to 
illustrate  the  concept  of  skewness  in  the  following  chapter,  and  to 
serve  as  a check  for  the  two-dimensional  models  far  from  the  sink. 

The  form  of  the  x-momentum  equation  solved  in  this  chapter  is 


Ie  „ _ 9xy* 

3x  3y 


(32) 


which  corresponds  to  laminar  one-dimensional  flow  with  no  body  forces. 

The  solutions  are  fairly  simple,  since  the  pressure  gradient  3p/3x  is 
constant  (and  is  represented  by  Px>.  The  equations  are  solved  first 
for  a two-layer  stratified  system,  then,  for  a linear  viscosity  gradient, 
and  finally,  for  a log  viscosity  gradient.  The  effect  of  a slip  boundary 
condition  is  then  studied. 


Two-Layer  Flow 


The  solution  to  the  stable,  two-layer  flow  of  immiscible 
fluids  problem  has  been  presented  by  Bird,  Stewart,  and  Lightfoot  (1960) 
as 


P * 

U - -5— 
T 8p_q 


UT  + WB 


+ - -■■  ■ (Tl)  - UYZ 

PT  + 


(33) 
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p ;.3 
u - -2— 

E 3pBq 


2,iB  UT  “ UP  7 

T + — 1 ( ?Y  ) - 4Y 

WT  + yT  + ya 


(34) 


where 


UT  = velocity  in  top  layer 
Ug  = velocity  in  bottom  layer 
y^  = viscosity  in  top  layer 
y„  55  viscosity  in  bottom  laver 

D 

Since  is  usually  not  known,  it  can  be  eliminated  by  ince- 
giating  the  velocity  profile  and  solving  for  in  terms  of  q,  Jt,  and  y. 


96(yfi  + yT)q  f 


PB  UT , 

PB2  + 14yByT  + yT2J 


(35) 


Substituting  back  into  the  velocity  equation  gives 


where 


12(y. 


uk  LPb  T Vv 


(36) 


WBPT 


+ 14V*Wt  + U 


B T mT 


Jpb»  for  bottom  layer 
K |y_,  for  top  layer 


A plot  of  this  solution  for  values  of  2(y  - y_)/(y_  + y ) = — 

B i B T — 

0 Jn  P 

which  corresponds  to  - — ^ for  this  problem  is  given  in  Figure  5. 

*o  y 

The  important  thing  to  note  is  that  as  the  viscosity  difference 
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between  layers  becomes  larger  more  fluid  is  withdrawn  from  the  top 
layer  and  the  flow  becomes  skewed.  For  the  case  of  Ap  « 0,  the  profile 
became  that  of  laminar  flow  In  a duct. 


Linear  Viscosity  Gradient 


For  the  case  of  a linear  viscosity  gradient 


JL-m  i + L.  IE  Y 

U 3v 
o 


3y 


(37) 


the  equation  of  motion  reduces  to 

.2., 


(38) 


where 


-p  “ Px*  /qp0 


P>0 


A « 2L  IE 
% 


p_  3y 


, 0<A<4 

with  no  slip  boundary  conditions 


U(±l/2)  - 0 
+1/2 

q‘  / 
-1/2 


U(Y)dY 


The  equation  can  be  solved  easiest  by  replacing  the  derivatives 
with  finite  differences.  The  equation  can  be  solved  by  estimating  a 
value  of  P,  solving  for  U(Y) , then  using  the  continuity  boundary  con- 


* 


1 


dition  to  correct  P.  This  can  be  repeated  until  P becomes  constant. 
An  analytical  solution  can  be  obtained  for  the  case  when  A - 0 (con- 
stant viscosity).  In  this  case  P - 12.  This  is  used  as  an  initial 
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estimate  of  P. 


The  finite  difference  equation  becomes 


N 

P - Pq/  5-  (U  AY) 
j - 1 3 

The  solution  is  plotted  for  three  values  of  A (.002,  2.,  3.98) 
to  show  the  effect  of  A on  the  velocity  profiles  in  Figure  6.  Figure  7 
shows  that  P will  remain  constant  for  small  values  of  A but  as  A becomes 
large,  P decreases  rapidly. 


Log  Viscosity  Gradient 


Some  runs  in  the  laboratory  portion  of  the  study  were  found  to 
have  viscosity  gradients  given  by 

V m Vi Q e”cly  (40) 

which  will  be  called  the  "log"  viscosity  gradient.  Substituting  for 
V and 


in 

ay 


-c,y 

~clyoe  1 


gives 


d2U  dU 
dY2  ' dY 


-p 


CY 


(41) 


(42) 


where 


C - c±l 

P 3 Px*3/qwo 
with  boundary  condition 
U(±l/2)  - 0 

The  solution  to  the  differential  equation  is 
ry 

u * ai  [a2  - e (a3  - 2Y) ] (43) 


where 


- -P/2C 

a2  = eC/2  [coth  (C/2)  -1] 

= coth  (C/2) 

Velocity  profiles  are  plotted  in  Figure  8 for  several  values  of  C. 
Since  P is  generally  not  known  (very  seldom  are  px  and  q both  known) , 
P was  determined  as  a function  of  C (C  / 0)  as 


P 


2C“ 


4 sinh  (C/2)  - c a. 


(44) 


A plot  of  this  function  is  given  in  Figure  9.  As  in  the  linear  case, 
P 12  as  C -*■  0. 


Slip  Along  Surface  Boundary 

A boundary  condition  which  may  be  of  some  practical  value  is 
that  in  which  the  flow  slips  along  a top  boundary  while  having  no  slip 
on  the  bottom.  This  type  of  boundary  condition  might  characterize 
flow  upstream  of  a weir. 


Fig.  8.  Velocity  profiles,  log  viscosity  gradient 
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The  x-momentum  equation  can  be  written 


dy  “ “X 

with  boundary  conditions 

Hf(+l/2)  - 0 = Tyx(+ 1/2 ) (46) 

which  is  equivalent  to  the  equation  solved  earlier  with  the  substitution 
of  the  dimensionless  Newton's  Law  of  Viscosity 

(47) 

U dY  (H,) 

o 

The  dimensionless  shear  stress  is  related  to  shear  stress  by 
*2 

T = T (48) 

qyQ  yx 

The  differential  equation  for  x can  be  solved  with  the  top  boundary  con- 
dition to  give 

x = -P/2(2Y  - 1)  (49) 

At  this  point,  Newton's  Law  of  Viscosity  can  be  substituted 


to  give 


_ du  = p ( i - 2y] 

% dY  Pl  2 J 


which  can  be  solved  with  p/pq  = e and  the  bottom  boundary  condition 
to  give 

PeCT  ?e'C/2 

U - ~~  (2YC  - C - 2)  + (1  + C)  (51) 

2CT  or 

A condition  for  continuity  can  be  used  to  derive  an  expression  for  P 


UdY  = 1 


CHAPTER  V 


SKEWNESS  OF  VELOCITY  PROFILES 

Definition 

As  was  shown  in  the  previous  chapter,  a variation  of  viscosity 
through  a fluid  tends  to  result  in  higher  velocities  in  regions  of 
lower  viscosity.  The  viscosity  tends  to  skew  the  velocity  distribution. 
The  concept  of  skewness  about  a point  refers  to  the  amount  of  a quantity 
that  lies  to  the  positive  or  negative  side  of  that  point. 

There  is  a need  in  this  study  for  an  unarbltrary  index  of  the 
effect  of  a viscosity  distribution  on  a velocity  profile.  As  was 
Implied  in  the  previous  paragraph,  the  concept  of  skewness  provides 
such  an  indicator.  The  most  conmonly  used  measure  of  skewness  is  the 
third  moment  of  the  distribution  about  a point.  The  third  moment  can 
be  defined  for  grouped  (discrete)  data  as 


where  is  the  point  about  which  the  moment  is  taken  and  f,  is  the 
° J 

frequency  (intensity)  of  whatever  is  distributed  about  Y . In  this 

o 

study  f represents  velocity.  Yq  is  taken  as  the  centroid  of  the 
velocity  profile  when  the  profile  is  symmetric.  This  corresponds  to 
the  elevation  Y ■ 0 throughout  this  study.  Skewness  has  dimensions  of 
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length  cubed.  To  avoid  problems  caused  by  scaling,  the  third  moment 
is  calculated  only  for  dimensionless  profiles  in  this  study. 

With  the  above  conventions,  the  formula  for  calculating  skew 


is  given  as 


i uj  v 

i-i  3 3 

N 

Z »j 

j - 1 3 


where 


Uj  - horizontal  velocity  at  elevation  Yj 
N = number  of  points  at  which  velocity  is  calculated 
When  the  velocity  profile  is  available  as  a continuous  function  U(Y), 


the  skew  can  be  defined 


UfY)Y3dY 


U(Y)dY 


Since  the  integral  in  the  denominator  represents  dimensionless  flow, 
which  will  always  be  one  since  U - ui/q,  the  definition  becomes 


skew  - J U(Y)YJdY 
-1/2 

Since  most  of  the  solutions  presented  in  this  report  are 
reached  using  numerical  methods,  the  skew  formula  using  discrete 
values  for  U is  used  most  often. 


Application  to  One-Dimensional  Velocity  Profiles 


The  velocity  profiles  developed  in  the  previous  chapter  all 
showed  Increasing  skew  as  the  viscosity  stratification  i/po(3p/8y)  term 
increased.  The  dimensionless  skewness  is  plotted  vs.  a parameter  that 
indicates  the  magnitude  of  the  viscosity  stratification  for  two-layer, 
linear  and  log  viscosity  gradients  in  Figures  11,  12,  and  13. 

The  figures  show  that  skewness  increases  with  the  appropriate 
viscosity  stratification  parameter.  There  exists  an  upper  bound  for 
the  parameter  in  both  the  two-layer  and  linear  case  since  the  change 
in  viscosity  cannot  exceed  twice  the  average  viscosity.  The  reader 
can  gain  an  appreciation  for  skewness  from  Figure  14  which  gives  the 
skewness  for  several  velocity  profiles. 

Application  to  Two-Dimensional  Velocity  Profiles 

In  two-dimensional  flow  to  a line  sink,  the  skew  increases  as 
one  moves  away  from  the  sink  where  skewness  is  zero.  It  is  therefore, 
necessary  to  specify  the  dimensionless  location  in  the  flume  at  which 
the  skewness  is  measured  in  the  two-dimensional  case. 

In  the  problems  which  require  upstream  boundary  conditions, 
these  conditions  will  also  play  a role  in  determining  skewness.  Most 
of  the  two-dimensional  problems  are  formulated  as  initial  value  prob- 
lems so  that  upstream  boundary  conditions  are  not  usually  required. 


CHAPTER  VI 


LINEAR  VISCOSITY  GRADIENT  NEGLECTING  BUOYANCY  FORCES 

Two  models  for  describing  viscosity  stratified  flow  are  pre- 
sented in  this  chapter.  The  first  model  is  a duct  model  which  is  used 
to  illustrate  some  concepts  for  two-dimensional  flow,  to  demonstrate 
the  importance  of  the  viscosity  stratification  term  in  the  x-momentum 
equation,  and  to  study  the  effect  of  including  inertia  terms  in  the 
model.  The  second  model  is  a tank  model  which  is  used  to  investigate 
the  effects  of  a fixed  back  wall  on  the  flow. 

Duct  Flow 

This  model  uses  the  boundary  layer  type  assumption  and  the 
assumption  that  buoyancy  forces  and  vertical  velocity  are  negligible 
to  reduce  equations  13  and  5 to 

iE  u_  ( 32U  . 32U|  + t_  30  3U  (59) 

"o  S*  “o  w n2J  "o  * SY  ' \ 9Y 
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Two-dimensionless  groups  are  apparent  in  the  above  equation, 
the  Reynolds  number  (qpo/pQ)  which  Indicates  the  magnitude  of  the 
Inertial  forces  and  the  viscosity  stratification  number  <*v/9y) 

which  indicates  the  importance  of  viscosity  stratification.  To  be 

consistent  with  its  uses  later  in  the  report,  this  number  is  defined  as 

51 
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A » i_  hL 
2 " u 3y 


The  equation  becomes 
P - i+|y 


. 32U  \ . A 3U 

3Y2 J + 2 3Y 


where 

p.a!|£ 
y 3x 
o 


with  boundary  conditions 
U(X,  ± 1/2)  » 0 

U(0,  Y)  - 6(0)  = Dirac  delta  function 


U(X  Y)  - 1 
max. 


1/2 


/ 


U(X,  Y)dY  - 1 


-1/2 


(61) 


(62) 


(63) 

(64) 

(65) 

(66) 


The  flow  starts  as  uniform  flow  upstream  and  flows  out  through 
a sink  at  Y * 0.  While  the  upstream  boundary  conditions  and  the  neglect 
of  buoyancy  forces  are  somewhat  unrealistic,  this  makes  the  model  easy 
to  use. 

Note  that  there  is  no  skewness  at  either  the  upstream  or  down- 
stream boundary.  With  the  lack  of  skew  upstream,  the  skewness  appears 
to  die  off  near  the  upstream  boundary.  This  is  an  anomoly  of  this 
model  which  is  dealt  with  later. 


Creeping  Flow  Results 

The  model  was  first  solved  with  the  inertial  terms  omitted 
using  the  techniques  described  earlier;  that  is,  solving  equation  59 


V a rxr 
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for  U(Y)  then  using  the  continuity  boundary  condition  to  get  a new  value 
for  P. 

The  first  effect  studied  was  that  of  omitting  the  A/2  3U/3Y  term 

and  allowing  the  varying  viscosity  to  be  accounted  for  in  the  (1  + AY/2) 

2 _ 

V U term.  The  results  are  shown  for  x/x  = X = 1/2  in  Figure  15. 

max 

Note  that  the  skew  is  almost  twice  as  large  for  the  case  in  which  the 

viscosity  stratification  is  accounted  for  using  both  terms.  The 

A/2  3U/3Y  term  will  be  retained  for  the  remainder  of  the  study. 

The  relationship  of  skew  vs.  distance  from  sink  is  shown  for 

two  values  of  A in  Figure  16.  The  skewness  increases  with  A much  as 

in  the  one- dimensional  cases  except  that  very  near  the  sink  the  profile 

becomes  symmetric.  This  accounts  for  the  decreased  skew  near  the  sink. 

The  model  was  also  run  for  a variety  of  length  to  depth  aspect 

ratios.  The  profile  for  different  values  of  X were  virtually  identical 

as  long  as  they  were  calculated  away  from  the  upstream  boundary.  For 

the  remaining  model  runs  in  this  section  the  aspect  ratio  x /£  will 

max 

be  set  to  10,  since  this  is  consistent  with  the  assumption  V«U  and 
yet  should  not  require  numerous  or  large  grid  cells. 

The  velocity  profile  development  from  the  uniform  flow  upstream 
to  the  narrow  withdrawal  layer  at  the  sink  is  shown  in  Figure  17.  Note 
that  once  the  flow  is  away  from  the  symmetric  upstream  and  downstream 
boundary  condition,  skewness  arises.  In  Figure  17  the  upstream  boundary 
is  at  x * 10. 

Effect  of  Reynolds  Number 

The  Reynolds  number  was  Increased  in  the  model  to  account  for 


k 

i 
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15,  Effect  of  3u/3y  term 


y ** 


Fig.  17.  Development  of  velocity  profile-neglecting  buoyancy 
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inertial  forces.  It  was  found  that  for  Reynolds  numbers  less  than  100, 
inertial  forces  had  no  effect  on  the  flow.  It  was  only  as  the  Reynolds 
number  approached  1000  that  the  profiles  began  to  change  dramatically. 
The  velocity  profile  became  blunter  and  the  skewness  decreased  dramati- 
cally. Because  of  this  result,  it  was  concluded  that  the  Reynolds 
numbers  to  be  investigated  in  this  study  would  be  0(100)  or  less.  In 
general,  the  laboratory  runs  would  be  for  Reynolds  number  between  10 
and  0.1. 


Tank  Boundary  Conditions 

The  equations  were  then  solved  with  boundary  conditions  corre- 
sponding to  that  of  a closed  end  tank  with  fixed  back  wall.  In  this 
case  the  boundary  conditions  were 


U(0,Y)  = 6(0) 

(67) 

U(X,  ± 1/2)  =*  0 

(68) 

U(X  Y)  - 0 

max. 

(69) 

q 


x 


U(X,Y)dX 


1 


x 

x 

max 


(70) 


The  solution  in  terms  of  streamlines  is  shown  in  Figure  18. 

The  growth  of  the  velocity  profile  is  shown  in  Figure  19.  The  profiles 
show  more  skew  in  the  tank  case  since  the  artificial  upstream  boundary 
condition  of  uniform  flow  was  not  imposed. 

This  Indicates  that  there  should  be  skew  even  when  the  flow 
was  virtually  zero  at  the  back  wall.  The  importance  of  the  back  wall 
was  investigated  at  this  point  in  the  study  because  it  is  impossible  to 


'V  «f 


reamlines  for  tank  model 


Development  of  velocity  profiio  k model 
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I 

I 

I 

have  a steady-state  model,  which  accounts  for  buoyancy  and  has  a fixed 
back  wall  since  a fixed  back  wall  implies  that  flow  is  caused  by  the 

► 

u falling  water  level  which  In  turn  means  that  the  density  field  is 

changing  with  time. 


I 


I 


CHAPTER  VII 


LINEAR  VISCOSITY  GRADIENT  WITH  BUOYANCY  FORCES 


Buoyancy  forces  have  not  been  accounted  for  in  the  velocity  pro- 
files determined  in  Chapters  IV  and  VI.  The  equations  developed  in 
Chapter  III  are  solved  below  and  the  skewness  is  given  for  an  array 
of  viscosity  and  density  gradients.  Criteria  for  existence  of  skewness 
are  then  derived. 

Model  Development 

2 2 

Equation  26  can  be  rearranged  to  give  (for  3 y/3y  = 0) 


where 


4 3 

+ a i_l 

yo  3Y4  3Y3 


A 


M in 

V 3Y 
o 


0 


B 


4 

8*  epo 

qPo 


(71) 


The  first  term  in  the  above  equation  represents  viscous  forces, 
the  second  viscosity  stratification  and  the  third  buoyancy.  The  dimen- 
sionless group  A represents  the  ratio  of  viscosity  stratification  to 
viscous  forces.  The  group  B can  be  written 


Reynolds  Number 

(Densimetric  Froude  Number) 


2 


(72) 
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which  is  actually  a balance  between  buoyancy  and  viscous  forces. 

To  solve  the  above  equation,  four  boundary  conditions  on  Y and 
one  on  X must  be  specified.  The  boundary  condition  on  X is  that  all 
flow  passes  through  the  origin 
-1/2  , Y>0 

Y(0,Y)  = ] 0 , Y-0  (73) 

1/2  , Y<0 


Two  conditions  on  Y are  used  to  insure  continuity 
Y(X,  ± 1/2)  - + 1/2 


(74) 


The  other  two  boundary  conditions  on  Y will  be  used  to  indicate  whether 
or  not  the  flow  is  slipping  along  the  top  and  bottom  of  the  duct.  They 
are 


■Mr  <x»  * 1/2)  ■ 0 » Blip 

3Y* 

^ u/ 

jY  (x»  1 1/2)  “0  , no  slip 


(75) 


The  above  equations  are  solved  numerically  as  described  in  Chapter  III. 

Figures  20,  21,  and  22  give  the  solution  to  the  slip  problem 
for  different  values  of  B.  Note  that  as  B increases,  the  withdrawal 
layer  becomes  smaller  as  buoyancy  forces  tend  to  inhibit  viscous 
momentum  transport.  Thin  can  be  deduced  from  the  equations  by  noting 
that  as  B becomes  large,  the  solution  reduces  to 


(76) 


which  corresponds  to  an  infinitely  thin  withdrawal  layer  with  infinite 
velocity. 
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As  B becomes  small,  the  solutions  approach  those  of  the  previous 
chapter.  As  B becomes  small  for  A = 0,  the  solution  corresponds  to 
uniform  flow  upstream  for  slip  boundary  conditions  and  to  a parabolic 
velocity  profile  for  no-slip  conditions. 

Equation  71  was  solved  for  slip  and  no-slip  boundary  conditions. 
Since  the  effects  of  viscosity  stratification  are  much  more  pronounced 
for  the  no-s] ip  boundary  conditions,  those  conditions  will  be  employed 
for  the  duration  of  this  chapter. 

Development  of  Velocity  Profiles 

The  solution  to  the  equations  of  this  chapter  in  terms  of 

3 

velocity  profiles  in  a duct  for  A ■ 4 and  B * 10  is  given  in  Figure 
23.  The  solution  can  be  divided  into  three  regions  in  the  x-direction. 

Near  the  sink,  buoyancy  forces  predominate  and  the  thickness 
of  the  forward  flowing  withdrawal  layer  is  much  less  than  that  of  the 
duct.  This  is  followed  by  a transition  region  in  which  the  buoyancy 
and  viacous  forces  are  of  similar  magnitude.  In  this  region  the 
velocity  profile  begins  to  drag  along  the  boundaries. 

Far  from  the  sink  the  viscous  forces  predominate.  The  stream- 
lines become  parallel,  the  pressure  gradient  becomes  hydrostatic,  and 
the  velocity  profiles  approach  those  given  in  the  one-dimensional 
solutions. 


Skewness  of  Velocity  Profiles 

The  skewness  of  the  velocity  profile  is  somewhat  irrelevant  in 
the  buoyancy  region  since  it  represents  the  importance  of  viscosity 


68 


stratification  which  is  negligible.  Furthermore,  the  existence  of  a 
backflow  zone  can  result  in  meaningless  values  for  skew. 

Equation  71  was  solved  for  a large  array  of  values  of  A and  B. 
The  skewness  was  found  to  be  a function  of  A,  B,  and  the  dimensionless 
distance  from  the  sink  (X). 

The  skewness  of  the  profiles  for  A = 2 and  A = 4 are  given  in 
Figures  24  and  25,  respectively.  Note  that  for  small  values  of  X and 
large  values  of  B the  buoyancy  forces  tend  to  mask  the  skew,  even  though 
A is  sufficiently  large  to  cause  significant  skew.  As  X becomes  large 
or  B small,  the  skewness  approaches  that  calculated  for  one-dimensional 
flows  which  correspond  to  the  case  B = 0. 

Figure  26  summarizes  the  results  of  these  solutions.  It  was 
decided  that  skew  became  significant  when  the  flow  above  the  line  Y = 0 
was  20  percent  more  than  that  below  the  line.  This  corresponded  to 
skew  of  0.002.  Figure  27  shows  this  velocity  profile,  which  occurs 
when  A = 1.2. 

For  A<1.2,  the  skew  is  negligible.  For  A>1.2  the  skew  becomes 

significant  far  from  the  sink  but  is  masked  by  buoyancy  forces  near  the 

sink.  The  transition  between  these  two  zones  depends  on  the  distance 

3 

from  the  sink.  For  X ■ 0(1),  B need  only  be  0(10  ) to  mask  skewness, 
while  for  X = 0(10),  B must  be  0(10^*). 

Implications 

It  is  clear  from  the  preceeding  discussion  that  viscosity 
stratification  can  occur  only  for  special  combinations  of  A,  B,  and  X. 

In  order  to  derive  a simple  formula  to  determine  if  skew  is  greater  than 
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0.002,  A,  B,  and  X were  plotted  to  gain  insight  into  their  interrela- 
tionship for  the  critical  value  of  skew. 

For  a constant  value  of  A,  both  B and  X are  plotted  as  straight 
lines  on  log-log  paper.  A multiple  regression  program  was  used  to 
determine  the  equation  for  these  lines  and  it  was  found  that 

B = 0.572(A-1.2)1/2X3/4  (77) 

fit  the  data  with  index  of  determination  0.895. 

To  give  additional  insight  into  the  criteria  for  viscosity 
stratification,  the  definitions  of  A,  B,  and  X were  substituted  into 
equation  77  and  the  equation  was  solved  for  3p/3y.  The  criteria  for 
viscosity  stratification  to  occur  is 


The  above  relation  shows  that  viscosity  stratification  will 
occur  for  small  values  of  the  density  gradient,  and  small  ducts  or 
at  large  distances  upstream  of  the  sink,  as  is  shown  in  the  laboratory 
results  in  Chapter  XI. 
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CHAPTER  VIII 


LOG  VISCOSITY  GRADIENT  WITH  BUOYANCY  FORCES 


In  this  chapter,  the  motion  equations  are  solved  for  the  log 
viscosity  gradient.  From  the  solution,  criteria  for  the  existence  of 
viscosity  stratification  are  developed. 


Model  Development 


As  described  earlier,  the  viscosity  gradient  for  the  log  vis- 
cosity profile  can  be  given  by 

y = Po  e~cly  (7 

where  c^  is  the  "slope"  of  the  gradient  on  semi-log  graph  paper,  and 
P is  the  viscosity  at  y * 0.  Defining  the  dimensionless  number 
C = c^i,  the  dimensionless  viscosity  and  its  derivatives  become 


r - A lE 

or  C * « 

u 3y 


y/P0  - e 


JL  | U „ _Ce-CT 

3y 


2 2 

* * V _ r2  -CY 

“ C.  e 

Po  3y2 


Substitution  of  the  above  derivatives  into  equation  26  from 

CY 

Chapter  III  yields,  after  multiplying  through  by  e 

_ ,r  ifl  + r2  2-1  + r *CY  W n 

4 2C  3 + C 2 + Be  "ay  ~ 0 
3Y  3Y  3Y  dX 
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As  in  the  case  of  the  linear  viscosity  gradient,  the  first  term 
represents  viscous  forces,  the  next  two  viscosity  stratification  and  the 
final  term  buoyancy. 

Equation  83  was  solved  with  the  no-slip  boundary  conditions 


||  (X,  ± 1/2)  = 0 
Y(X,  ± 1/2)  = + 1/2 


-1/2  , Y>0 

Y(0,Y)  = ^ 0 , Y=0 

1/2  , Y<0 


(84) 

(85) 

(85) 


The  velocity  profiles  behaved  similarly  to  those  with  linear 
gradients  (i.e.,  buoyancy  forces  predominating  near  sink,  viscous 
forces  predominating  far  from  sink).  For  B = 0,  the  solution  approached 
that  obtained  for  the  one-dimensional  case.  This  is  shown  in  Figure  28. 

Skewness  of  Profiles 


The  skew  of  the  velocity  profiles  were  calculated  for  an  array 
of  values  for  B,  C,  and  X.  These  are  presented  in  Figures  29,  30,  and 
31.  The  skewness  approached  that  of  the  one-dimensional  solution  for 
large  values  of  X and  small  values  of  B. 

A plot  of  B vs.  C showing  zones  in  which  viscosity  stratifica- 
tion can  occur  is  given  in  Figure  32.  Combinations  of  C,  B,  and  X 
that  yield  values  of  critical  skew  (0.002)  were  related  by 

B - 555(C  - 0. 54)°' 30  X 1,03  (87) 


with  index  of  determination  0.988.  Approximating  the  exponent  on  X as 
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CX0I5 


Fig.  30.  Skew 
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1.0  and  substituting  for  the  definitions  of  B,  C,  and  X gives  the 
criteria  for  skewness  as 


> a 


+ 0.54 


(88) 


Significant  skewness  will  therefore  occur  for  large  distances  from  the 
sink  when  the  viscosity  gradient  is  large. 


Incipient  Skewness 

In  the  case  of  the  linear  viscosity  gradients,  the  flow  devel- 
oped from  being  buoyancy  dominated  to  viscous  dominated  over  a fairly 
short  transition  region.  In  the  case  of  a log  viscosity  gradient 
this  transition  zone  is  much  larger. 

The  existence  of  the  transition  zone  is  due  partly  to  the  fact 
that  the  fluid,  being  more  viscous  near  the  bottom,  "grabs  onto"  the 
bottom  boundary  much  sooner  than  it  grabs  onto  the  top  boundary.  There 
is  still  backflow  in  the  transition  zone  so  the  skewness  coefficient 
is  not  a good  indicator  of  the  effects  of  the  viscosity  stratification. 

This  skewing  of  the  profile  in  the  transition  zone  where  viscous 
forces  dominate  in  the  bottom  of  the  flume  and  buoyancy  forces  dominate 
in  the  top  will  be  defined  as  "incipient  skewness."  The  four  graphs 
in  Figure  33  show  how  the  flow  develops  from  buoyancy  dominated  with 
no  skewness,  to  incipient  skewness  and  finally  fully  developed  skewness. 

Figure  34  shows  schematically  how  the  zone  of  incipient  skewness 
relates  to  the  buoyancy  and  viscous  zones.  The  boundaries  between  the 
zones  are  not  intended  to  imply  that  this  is  their  actual  shape. 

Streamlines  are  plotted  in  Figure  35  to  illustrate  the 


Fig.  33.  Development  of  velocity  profile-log  gradient 
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Fig.  35.  Streamlines  for  Incipient  skewness 
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transition  from  incipient  skewness  to  fully  developed  skewness.  The 
eddy  at  the  surface  represents  the  backflow  that  occurs  in  the  buoyancy 
dominated  zone.  The  streamlines  become  horizontal  for  large  values  of 
X.  The  values  which  they  approach  are  shown  at  X = °°. 

Incipient  skew  begins  once  the  0.5  streamline  reaches  the 
bottom  boundary.  Using  this  as  an  indicator  of  incipient  skew.  Figure 
36  was  prepared.  It  gives  the  value  of  X at  which  this  occurs  for  an 
array  of  values  for  B and  C. 

The  multiple  regression  equation  which  fits  the  data  in  Figure 

36  is 

X = B • e_0,4C(3.7  x 10"5)  (89) 

Substituting  for  X,  B,  and  C gives 

x = 3.7  x 10  5 e 0,4C (£5gePQ/qpo)  (90) 

x can  be  interpreted  as  the  distance  upstream  of  the  sink  at  which 
separation  of  the  0.5  streamline  from  the  wall  occurs. 


CHAPTER  IX 


EFFECT  OF  UPSTREAM  BOUNDARY  CONDITIONS 
AND  INERTIA  ON  VELOCITY 


In  developing  the  mathematical  model  in  Chapter  III,  it  was 

assumed  chat  U y/y  could  be  replaced  by 
oo  o 

1 + AY/2  for  linear  gradient  / 


or 


-CY 

e for  log  gradient 

This  involved  two  assumptions:  U = 1 and  the  viscosity  field  is  fixed 
in  space.  The  effect  of  these  assumptions  and  of  the  existence  of  a 
fixed  back  wall  in  the  experiments  conducted  are  discussed  in  this 
chapter. 


Viscosity  Constant  on  Streamline 


While  in  a non-diffuse  model,  viscosity  is  strictly  a function 
of  the  stream  function  (V)  ; the  viscosity  has  been  modeled  thus  far  as 
a function  of  depth.  This  type  of  assumption  makes  mathematical 
modeling  of  the  system  much  easier,  but  may  make  the  results  inaccurate. 
The  correct  approach  to  modeling  the  viscosity  is  to  treat  the  vis- 
cosity as  constant  on  a streamline  (y  = f(T)  only). 

To  find  the  viscosity  at  a point,  it  is  necessary  to  know  the 
stream  function  at  that  point.  The  problem  there  is  to  set  up  a re- 
lationship between  Y and  y.  All  that  is  known  is  that  when  the  fluid 
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is  at  rest,  p/pq  =»  e CT.  (The  log  gradient  will  be  used  here  to  illus- 
trate the  implications  of  using  the  fixed  viscosity  field.) 

When  steady-state  is  reached,  there  should  be  some  point  at 
which  the  relationship  below  is  true.  This  is  taken  as  x >*  ® 

V/UQ  (x  - »)  - e~CT  (91) 

even  when  flow  is  occurring. 

It  is  necessary  then  to  derive  a relationship  between  ¥ and  Y 
at  X * ®.  Ideally , the  relationship  is  the  solution  for  one-dimensional 
flow.  Unfortunately,  those  equations  cannot  be  solved  for  Y ■ f(y). 

So  to  simplify  matters,  a uniform  upstream  velocity  profile  is  used 


which  implies 
¥ * -Y 

for  the  boundary  condition  ¥(0)  - 0. 

The  viscosity  at  any  point  in  the  fluid  can  be  given  as 


(92) 


(93) 


( CY 

e T , -5 <¥<0.5 

U/U0(Y) 

I eC(±0. 5)  , ¥<±0.5 


(94) 


Equation  26  with  this  new  expression  for  viscosity  becomes 


2c  + r2  + r .-CT  3¥ 

2C  , + C -jHe  a?  " 0 


3Y' 


3Y 


(95) 


This  equation  was  solved  numerically  and  the  solution  is 
plotted  in  Figure  37,  along  with  the  solution  from  the  previous  chapter. 
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Note  that  the  solutions  do  not  differ  by  much  even  though  this  solution 
is  for  the  incipient  skewness  zone  where  one  would  expect  the  profile 
to  be  most  sensitive. 


Non-Linear  Upstream  Velocity  Profile 


For  all  solutions  given  thus  far,  it  was  assumed  that  at  some 
point  very  far  from  the  sink,  that  the  velocity  was  uniformly  distri- 
buted. This  permitted  writing  the  stream  function  at  that  point  as 

4-  = -Y  (g 

With  this  assumption,  it  was  possible  to  define  simple  relationships 
between  density,  viscosity  and  stream  function,  namely 


3(PS  + P) 


iie.il 


since 


These  assumptions  were  removed  and  the  upstream  velocity  profile 
was  taken  from  the  one-dimensional  solution  for  the  log  gradient  as 


where  a^...a^  are  defined  in  Chapter  IV.  It  is  now  no  longer  possible 
to  write  simple  expressions  for  p and  p as  a function  of  ¥. 

The  equations  were  therefore  solved  by  first  setting  up  a table 
of  p and  p vs.  V upstream.  Then  after  each  iteration  a value  of  vis- 
cosity and  density  could  be  assigned  to  each  grid  cell  by  interpolating 
from  the  internal  table. 


Equation  95  was  solved  and  the  results  were  found  to  be  identi- 
cal in  virtually  all  cases.  This  could  have  been  anticipated  by  noting 
that  the  corrections  on  both  viscosity  and  density  were  in  the  B term  in 
the  equations.  From  earlier  solutions  it  was  clear  that  B had  to  be 
changed  a great  deal  to  significantly  affect  the  solutions.  The  correc- 
tions introduced  in  this  section  only  changed  B slightly  from  the  values 
obtained  from  the  cssumption  of  a linear  profile  far  upstream. 

The  principal  difference  between  solutions  developed  from  this 
method  and  those  from  methods  presented  earlier  is  the  cost  of  the  com- 
putation, which  increased  greatly  because  of  the  table  lookups.  This 
solution  method  was  not  investigated  any  further. 

Effect  of  Back  Wall  on  Density  Field 

In  a tank  with  a fixed  back  wall,  it  is  not  possible  to  maintain 
a constant  density  field  since  layers  of  the  fluid  near  the  sink  eleva- 
tion are  being  withdrawn  much  as  one  would  withdraw  a single  card  from 
a deck  of  playing  cards.  To  correctly  model  a tank,  an  unsteady-state 
model  is  necessary  to  account  for  the  changing  density  field.  The 
question  is  one  of  how  different  are  the  flows  in  an  infinite  duct  with 
a steady-state  model  from  those  in  a tank. 

Since  in  this  experiment  the  flow  is  such  that  the  surface  of 
the  tank  drops  extremely  slowly,  the  velocity  field  similarly  changes 
very  slowly.  The  principal  effect  is  that  in  the  tank  the  density  is 
not  constant  along  a streamline  as  it  is  in  the  duct  model. 

The  species  (or  density)  continuity  equation  for  unsteady, 
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non-diffuse  flow  is 
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While  3p/3t  can  be  known  everywhere  only  by  solving  the  complete  unsteady 
momentum  and  continuity  equations,  it  can  be  stated  that  at  the  top 
boundary 
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where  V*  is  the  velocity  at  which  the  surface  is  dropping, 
boundary 


At  the  bottom 
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The  effect  of  this  term  on  the  solution  can  be  seen  by  remem- 
bering that  the  species  continuity  equation  entered  the  mathematical 
model  by  replacing  3p/3x  by  Vep  /U  . The  substitution  now  becomes 
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Since  3p/3x  is  substituted  into  the  buoyancy  force  term,  it 
can  be  concluded  that  the  existence  of  the  back  wall  of  the  tank  reduces 
the  effect  of  the  buoyancy  forces.  The  magnitude  of  this  effect  depends 
on  the  length  to  depth  aspect  ratio  of  the  tank.  If  it  is  small,  the 
effect  of  the  falling  surface  can  be  significant.  The  experimental  data 
in  this  study  indicate  that  buoyancy  forces  were  somewhat  smaller  than 
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predicted.  This  may  be  explained  by  the  moderate  length  to  depth  aspect 
ratio  used. 


Effect  of  Inertial  Forces 


In  most  of  the  solutions  presented  thus  far,  inertial  forces 
have  been  neglected.  The  experimental  results  presented  in  Chapter  XI 
indicate  that  this  did  not  seriously  affect  the  correlation  of  results 
from  the  mathematical  and  experimental  model,  since  the  Reynolds  numbers 
were  small  in  the  experimental  model. 

The  question  remains  as  to  what  the  mathematical  model  would 
look  like  if  inertial  forces  were  included.  Such  a model  using  similar 
assumptions  to  those  presented  in  Chapter  III  can  be  written 
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(viscous)  (viscosity  stratification)  (buoyancy)  (inertia) 
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where  F2  = R/B  = q2/i^eg 

Near  the  sink  the  inertia  and  buoyancy  terms  tend  to  mask  viscous 
forces  but  as  X ■*  the  velocity  profiles  will  approach  those  presented 
for  one-dimensional  flow.  This  viscous  zone  most  likely  occurs  even 
farther  upstream  when  R becomes  larger.  When  buoyancy  forces  are  large 
(small  F) , they  not  only  can  mask  the  viscous  forces  but  the  inertial 
forces  as  well.  The  buoyancy  forces  may  also  inhibit  turbulence  which 
will  make  the  above  model  applicable  at  much  higher  Reynolds  numbers 
than  if  F were  large. 
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CHAPTER  X 


LABORATORY  PROCEDURES 

The  velocity  profiles  were  measured  In  experimental  flumes  by 
photographing  the  locations  of  dye  streaks  at  various  times.  The  de- 
sired gradients  of  viscosity  and  density  were  developed  by  placing 
layers  with  decreasing  sugar  or  CMC  (carboxymethylcellulose)  concentra- 
tion in  the  flume.  This  chapter  of  the  report  describes  the  equipment 
and  procedures  used  in  conducting  the  experiments. 

Experimental  Flumes 

The  physical  model  runs  were  conducted  in  three  different 
flumes.  The  three  flumes  are  referred  to  in  the  report  as  the  big, 
narrow,  and  wide  flumes.  They  are  pictured  in  Figures  38,  39,  and  40, 
and  their  dimensions  are  given  in  Table  1.  The  length  is  the  distance 
from  the  sink  to  the  back  wall  of  the  flume;  the  height  is  the  distance 
from  the  bottom  to  the  greatest  depth  used;  and  the  width  is  the  distance 
between  the  sides.  The  camera  location  is  the  distance  from  the  sink  to 
the  point  at  which  the  camera  recording  the  position  of  the  dye  streaks 
was  usually  located.  This  procedure  is  described  later. 

The  big  flume  was  used  first.  When  the  density  gradients  used 
were  large,  the  results  of  the  experimental  runs  were  satisfactory. 

When  the  density  gradients  became  small,  convective  mixing  currents 
in  the  tank  became  significant  since  temperature  was  not  controlled  in 
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Fig.  40.  Wide  flume 
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TABLE  1 

CHARACTERISTICS  OF  FLUMES 


Flume 

Length  (cm) 

Height  (cm) 

Width  (cm) 

Camera 

Location  (cm) 

Big 

732 

60 

30 

73,  222,  317 

Narrow 

78 

58 

3 

32 

Wide 

67 

12 

30 

31 

the  building  housing  the  flume.  (It  would  change  by  as  much  as  15°C 
during  the  course  of  the  day.)  These  secondary  currents  often  inter- 
fered with  the  experiments. 

Furthermore,  large  concentration  gradients  could  not  be  used 
because  of  the  costs  and  difficulty  of  handling  large  quantities  of  the 
stratifying  fluid.  For  example,  it  would  require  one-half  tonne  of 
sugar  to  achieve  some  of  the  desired  gradients. 

The  narrow  flume  was  then  constructed  so  that  the  experiments 
could  be  run  in  a controlled  environment  and  would  not  require  large 
quantities  of  the  model  fluids.  Most  of  the  successful  runs  were  con- 
ducted in  this  flume.  Some  question  arose  as  to  possible  sidewall 
effects  in  such  a narrow  flume.  To  avoid  the  effects  of  the  side- 
wall,  the  wide  flume  was  constructed.  This  flume  had  a depth  to  width 
aspect  ratio  of  1:3. 

In  this  study  the  depth  to  width  ratio  was  varied  from  1:3  to 
20:1.  This  was  a much  wider  range  than  previous  Investigators,  as 
Koh  (1964)  used  3:1,  Imberger  (1970)  used  4:1,  and  Kao  et  al.  (1972) 


TT  * 


99 


used  2:1. 

The  flow  was  regulated  in  the  large  flume  by  a gate  in  the  down- 
stream wall  in  the  big  flume.  The  opening  in  the  gate  was  controlled  by 
a set  of  screws  at  the  top  of  the  gate  (see  Figure  38) . In  the  narrow 
and  wide  flumes,  the  flow  was  regulated  by  a clamp  on  the  tubing  flowing 
out  of  the  sink  (see  Figures  39  and  40). 

Measuring  Velocity  Profiles 

The  velocity  profiles  were  measured  by  determining  the  displace- 
ment of  dye  streaks  during  known  time  intervals.  Once  the  flow  had 
reached  steady-state,  small  pellets  of  erio-glaucine  dye  were  dropped 
into  the  flume.  These  were  photographed  at  several  intervals.  A clock 
was  placed  in  the  photograph  to  record  the  time  at  which  each  was  taken. 

The  pictures,  which  were  taken  as  35-mm  slides,  were  read  using 
a microfiche  reader.  The  location  of  the  dye  streak  and  time  of  the 
slide  were  recorded  and  input  into  a computer  program,  VELPRO.  VELPRO 
corrected  the  profiles  for  parallax  error,  calculated  the  velocity 
profile  and  integrated  the  profile  to  obtain  the  flow  rate.  The 
formulas  for  parallax  correction  and  a description  of  the  VELPRO  program 
are  given  in  Appendix  B. 

Figure  41  is  a schematic  diagram  of  a typical  dye  streak  during  a 
run.  Figures  42,  43,  and  44  show  actual  dye  streaks  in  the  big,  narrow, 
and  wide  flumes,  respectively.  Only  one  camera  was  used  for  the  narrow 
and  wide  flumes.  Three  cameras  were  used  in  the  big  flume  since  the  en- 
tire flume  could  not  fit  into  one  shot.  The  horizontal  distance  from  the 
sink  to  the  point  at  which  the  cameras  were  perpendicular  is  listed  for 
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Fig.  43.  Typical  photo,  narrow  flume 
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each  flume  in  Table  1 and  is  shown  for  the  large  flume  in  Figure  45. 

Model  Fluids 


The  desired  density  and  viscosity  gradients  were  achieved  in  the 
flume  by  filling  the  flume  with  layers  of  sugar  or  CMC  solutions.  While 
similar  viscosity  gradients  could  be  achieved  with  either  fluid,  the 
sugar  created  a much  larger  density  gradient  than  CMC  for  the  same 
viscosity  gradient. 

For  each  run  it  was  necessary  to  determine  the  viscosity  and 
density  at  every  point  in  the  fluid.  Since  it  was  impossible  to  make 
direct  measurements,  it  was  necessary  to  measure  some  property  of  the 
fluid  which  could  be  easily  tested  and  to  correlate  this  property  with 
the  viscosity  and  density.  For  sugar  this  property  was  refractance; 
for  CMC  it  was  conductivity.  In  Appendix  C,  these  correlations  are 
developed. 

Since  CMC  is  a generic  name,  there  are  many  types  of  CMC's 
available.  The  four  types  used  in  this  study  are  D-70  (detergent 
grade  from  the  Buckeye  Cellulose  Corp.),  7M,  7H,  and  7H4  (technical 
grade  from  Hercules  Incorporated) . The  D-70  CMC  contains  a great  deal 
of  suspended  cellulose  which  made  it  very  cloudy  and  tend  to  settle. 

D-70  CMC  was  used  only  in  a few  runs.  The  7M,  7H,  and  7H4  CMC's 
dissolved  completely  and  were  used  to  a greater  extent.  The  7M  CMC 
produced  moderate  viscosities  while  the  7H  and  7H4  produced  much  higher 
viscosities  at  the  same  concentration. 

In  high  concentrations,  CMC  behaves  as  a non-Newtonian  fluid. 
Since  most  of  the  theoretical  development  in  this  study  is  based  on 
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the  assumption  of  non-Newtonian  behavior,  it  was  important  that  the  CMC 
be  used  only  in  concentrations  and  at  shear  rates  for  which  it  is 
Newtonian.  In  Appendix  C,  it  is  shown  that  CMC  in  concentrations  less 
than  5 g/i  behaves  as  a Newtonian  fluid  in  a wide  range  of  shear  rates. 

Measuring  Concentration  Gradients 

Concentration  gradients  were  measured  after  each  run  so  that  the 

viscosity  and  density  gradients  could  be  determined.  Since  the  refrac- 

3 

tance  test  requires  only  1 cm  of  the  sample,  it  was  also  conducted 

before  each  sugar  run.  Since  the  conductivity  test  requires  approxi- 
3 

mately  100  cm  of  sample,  it  was  not  conducted  before  each  run  for  fear 
that  it  might  disturb  the  gradients. 

The  sugar  samples  were  taken  using  an  L-shaped  glass  tube  (see 
Figure  39) . The  CMC  samples  were  taken  from  a sampling  port  located 
at  the  upstream  end  of  the  narrow  and  wide  tanks. 

Flume  Filling  Procedure 

The  required  concentration  gradient  for  a run  was  determined 
from  the  desired  viscosity  and  density  gradient.  Then  the  numbers  of 
layers  to  be  input  to  the  flume  (varying  from  five  to  thirteen,  depend- 
ing on  the  flume)  were  selected.  Once  the  concentration  gradient  and 
number  of  l4yers  were  chosen,  the  concentration  in  each  layer  could  be 
calculated  from 


where 


n “ number  of  layers 
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3s/9y  = concentration  gradient 
i,  = depth  of  flume 
j = concentration  in  jth  layer 

The  flume  was  filled  starting  with  the  bottom  layer  using  the 
apparatus  shown  in  Figure  46  for  the  wide  and  narrow  flumes.  The  fluid 
was  poured  into  the  top  container.  The  flow  was  then  regulated  by  the 
elevation  of  the  top  container  and  a stopcock.  The  fluid  flowed  through 
the  tubing  to  the  plate  located  at  the  desired  elevation  of  the  layer. 

A sponge  was  located  on  the  plate  to  dissipate  any  momentum  the  fluid 
may  have  had  to  prevent  mixing  of  the  layers. 

In  the  case  of  the  big  flume,  the  fluid  was  pumped  in  through 
a hose  onto  a much  larger  plate.  Very  little  mixing  occurred  when  the 
narrow  and  wide  flumes  were  filled.  Some  mixing  did  occur  in  the  big 
flume  so  the  actual  concentration  gradients  were  somewhat  smaller  than 
the  desired. 

In  later  runs  the  layers  were  alternately  colored  orange  and 
clear.  This  was  done  to  better  track  motion  of  the  layers  and  to 
observe  when  the  layers  began  diffusing  into  one  another. 

Once  the  flume  was  filled  (which  took  approximately  one  day), 
from  one  to  three  days  were  allowed  to  elapse  to  let  the  interfaces 
between  the  layers  diffuse  into  a linear  gradient. 

Run  Procedure 

Before  the  run  the  concentration  gradient  was  checked  and  a dye 
particle  was  dropped  and  observed  to  determine  if  there  was  significant 
background  motion.  The  clock  and  run  number  were  set  up  in  front  of  the 


Fig.  46.  Flume  filling  technique 
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flume  and  the  camera  and  lighting  were  checked. 

The  flow  was  initiated  by  opening  the  clamp  or  gate,  and  the  rate 
was  checked  using  a cup  and  stopwatch.  Once  the  rate  was  set,  it  was 
usually  not  adjusted  as  transient  currents  might  have  interfered  with 
results. 


The  flow  was  allowed  to  come  to  steady-state, 
given  by  Koh  (1966b)  as 

t 

V^6s 


This  time  was 


(107) 


where 

6 = thickness  of  withdrawal  zone,  cm 

s * 

x = distance  from  sink,  cm 

For  the  big  flume,  this  time  was  on  the  order  of  5 minutes.  For  the 
narrow  and  wide  flumes  it  was  on  the  order  of  1 minute.  Usually,  approx- 
imately twice  this  time  was  allowed  before  velocities  were  measured. 

Once  steady-state  was  achieved,  the  dye  pellets  were  dropped. 
While  several  may  have  been  dropped  to  illustrate  withdrawal  layer 
growth,  one  was  always  dropped  at  the  location  labeled  "camera  location" 
in  Table  1.  The  dye  was  dropped  as  close  to  the  flume  center  line  as 
possible  to  minimize  any  side-wall  effects.  The  dye  streak  was  photo- 
graphed at  least  four  times  as  it  deformed. 

Once  the  run  was  completed,  the  flow  rate  might  have  been  ad- 
justed and  the  run  repeated.  In  the  narrow  flume,  only  three  runs  could 
be  attempted  until  the  surface  elevation  dropped  significantly.  In  the 
other  flumes,  this  effect  did  not  occur,  but  usually  only  three  runs 
were  made  with  any  flume. 


CHAPTER  XI 


EXPERIMENTAL  RESULTS 

The  results  of  experimental  runs  are  presented  in  this  chapter. 
First  the  parameters  for  the  runs  are  presented.  The  mathematical 
model  results  are  compared  with  the  experimental  results  on  a run-by-run 
basis.  The  overall  results  of  the  experimental  runs  are  summarized  in 
terms  of  the  dimensionless  numbers  B and  C. 

Data  Summary 

Table  2 summarizes  the  values  for  the  most  important  parameters 
for  each  run.  Each  column  is  explained  below.  Instruments  used  to 
measure  the  parameters  are  discussed  in  more  detail  in  Appendix  C. 

The  run  number  is  the  ordinal  number  of  the  run.  Runs  in  the 
big  flume  have  only  one  number  since  every  time  a parameter  was  changed , 
the  run  number  was  changed.  Runs  in  the  other  two  flumes  take  the  form 
M.N  where  M is  the  run  number  assigned  when  the  flume  is  filled, 
while  N is  the  number  of  the  flow  rate  used.  For  example,  run  50.1  is 
run  50  with  flow  rate  1 while  50.3  is  run  50  with  flow  rate  3. 

Note  that  many  run  numbers  do  not  appear  in  the  table.  Runs 
20-29  were  made  using  CMC  in  the  small  flume.  There  were  problems 
with  turbulence  in  most  of  these  runs  which  rendered  them  useless.  Some 
of  the  runs  with  numbers  in  the  30 's  were  attempted  while  additional 
model  fluid  was  injected  in  the  sampling  ports  at  the  upstream  end  of 
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the  narrow  flume  to  maintain  steady-state  conditions.  It  was  found 
that  inertia  from  these  Inflows  caused  turbulence  so  this  method  was 
abandoned . 

The  "flume"  column  contains  the  abbreviated  flume  name:  B-big 
N-narrow,  and  W-wide,  as  defined  earlier.  The  "model  fluid"  column 
designates  whether  the  model  fluid  was  sugar  (S)  or  one  of  the  CMC's 


(D-70,  7H,  7M  or  7H4) . 

2 

The  flow  at  the  sink  is  given  in  the  "q"  column  in  cm  /sec. 
This  value  was  determined  from  the  output  from  the  VELPRO  program  and 
the  fact  that  q increased  linearly  from  zero  at  the  back  wall.  This 
is  given  by  the  formula 


where 


qx  is  flow  at  position  x 

x is  length  of  flume 
max 

The  relationship  is  plotted  in  Figure  47  for  three  runs  which  had 
velocity  profiles  measured  at  several  locations. 

Po  and  are  the  density  and  viscosity  at  the  sink  elevation, 
e was  determined  by  plotting  density  vs.  depth  and  fitting  a straight 
line.  The  slope  of  the  line  was  Sp^/Sy  and  e « -l/po  3Pg/3y» 

l is  the  distance  from  the  bottom  of  the  fluid  to  the  surface 
of  the  model  fluid.  The  sink  was  always  at  approximately  £/ 2 from  each 
boundary.  Note  that  in  some  runs  a false  bottom  was  developed  by 
filling  the  bottom  of  the  flume  with  a very  dense  sugar  solution.  The 
CMC  or  lighter  sugar  solutions  were  laid  on  top  as  shown  in  Figure  48. 


q (cm/  sec; 


Fig.  47.  Flow  rate  as  function  of  distance  from  sink 
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This  allowed  l to  be  set  at  a much  smaller  value  than  the  total  tank 
depth. 


The  parameter  C was  defined  earlier  as  where  is  the  slope 

of  a semi-log  plot  of  viscosity  vs.  elevation.  The  parameter  B is 

the  ratio  of  the  Reynolds  number  (R  = qp  /y  ) and  the  densimetric  Froude 

o o 

number  (F  = q/i.  ~\fe g)  squared.  That  is 


4 

g l ep( 
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Note  that  the  B and  F are  written  in  scientific  notation  where 
4.2E+8  is  4.2xl08  and  2.3E-4  is  2.3xl0-4. 

The  experimental  runs  covered  a fairly  wide  range  of  the  im- 
portant parameters 

1.1  x 103  < B < 2.4  x 1011 


0.002  < C < 6.6 


0.054  £ R £ 156 

1.5  x 10~6  £ F £ 2.7  x 10"2 

All  the  runs  are  summarized  in  Figure  49  with  respect  to  the 
values  of  B and  C.  As  was  predicted  by  the  mathematical  model,  runs 
with  large  values  of  B do  not  show  skew  because  of  the  masking  effect 
of  the  large  buoyancy  forces.  As  B decreases  skewness  becomes  evident 
and  finally  predominates  for  sufficiently  low  B. 


Velocity  Profile  Results 

Figures  50-54  contain  graphs  of  the  observed  (experimental 
model)  velocity  profiles  plotted  with  the  corresponding  predicted 
(mathematical  model)  profiles  from  the  numerical  solution  of  equation 
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95.  The  actual  velocity  profiles  have  been  smoothed  In  some  cases 
where  there  was  turbulence  or  a secondary  surface  current.  For  example. 
Figure  42  shows  a forward  flowing  current  at  the  surface  (probably  due 
to  wind  stress).  This  current  is  not  shown  in  the  graph  of  run  11. 

Velocity  profiles  are  given  in  Figures  50-54  for  16  different 
runs.  They  are  grouped  according  to  the  size  of  the  flume  and  the 
magnitude  of  B.  The  values  of  B,  C,  and  X for  the  run  are  given  on  each 
figure. 

Figures  50  and  51a  show  runs  in  the  big  flume  with  large  values 
of  B (runs  8,  11,  12,  13,  14).  Note  that  there  is  a small  withdrawal 
layer  with  some  backflow.  Outside  the  withdrawal  layer  there  is  no 
horizontal  velocity. 

The  next  two  figures  (51b  and  51c)  are  for  runs  36  and  41  with 
similar  values  of  B in  the  narrow  flume.  They  show  essentially  the 
same  profiles  even  though  the  scales  of  the  experiments  were  smaller. 

Run  53  was  conducted  in  the  wide  flume  and  shows  the  same  profile 
shape  (Figure  51d) . 

The  preceding  runs  were  all  for  small  Reynolds  numbers  (0.58 
to  23).  To  test  the  applicability  of  the  model  at  higher  Reynolds 
numbers,  runs  54  and  56  were  conducted  at  Reynolds  numbers  of  85  and 
156.  The  model  still  predicts  the  profile  shape  as  shown  in  Figures 
52a  and  52b. 

From  the  preceding  profiles  it  is  c.  . ir  that  the  model  predicts 
velocity  profiles  well  for  runs  with  large  values  of  B.  The  model  was 
also  tested  for  smaller  values  of  B where  viscosity  stratification 
becomes  significant. 
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Runs  42  and  43,  in  Figures  52c  and  52d,  illustrate  the  occur- 
rence of  incipient  skewness  in  the  narrow  flume.  The  bottom  of  the  pro- 
file looks  like  the  one-dimensional  velocity  profiles  while  the  top 
resembles  the  buoyancy  dominated  profiles.  The  same  type  of  profiles 
exist  in  the  wide  flume  for  runs  55  and  57  (Figures  53a  and  53b). 

There  were  some  problems  with  run  57  since  in  the  mathematical 
model  prediction  the  profile  has  not  yet  begun  to  skew  while  in  the 
physical  model  significant  skewing  has  occurred.  While  this  is  more 
dramatic  in  this  run  than  some  others,  it  occurred  in  some  runs  with 
small  values  of  B.  That  is,  the  forward  flowing  layer  reaches  the 
boundaries  for  lower  values  of  X in  the  physical  model.  This  can  be  ex- 
plained by:  (1)  the  non-linear  inertial  terms  have  been  neglected  near 
the  sink.  Imberger  has  shown  that  the  predicted  withdrawal  layer  is 
larger  when  these  terms  are  included  and  would  therefore  match  the  ex- 
perimental data  better,  and  (2)  in  applying  a duct  model  to  experimental 
data  from  a tank  with  a fixed  back  wall,  the  horizontal  accelerations 
are  not  accounted  for  in  the  mathematical  model.  This  is  explained  in 
greater  detail  in  Chapter  IX. 

Nevertheless,  the  mathematical  model  does  account  for  the 
transition  from  buoyancy  dominated  to  viscosity  dominated  flow  very 
well.  This  is  illustrated  for  run  55  (Figure  53a)  in  that  the  experi- 
mental profiles  are  matched  with  predicted  profiles  both  near  the  sink 
and  away  from  the  sink  for  the  same  run. 

Fully  developed  skewness  was  observed  only  in  two  experimental 
runs  (48  and  50).  In  these  cases  the  experimental  profile  can  actually 
be  matched  with  the  one-dimensional  profile.  Runs  48  and  50  (Figures 
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53c  and  53d)  show  no  backflow  at  the  surface  while  the  peaks  in  the 
observed  and  predicted  profiles  occur  at  the  same  elevation. 

The  peaks  of  the  profiles  were  considerably  higher  in  the 
experimental  profile.  It  was  postulated  that  this  was  due  to  the 
fact  that  the  CMC  did  not  have  time  to  form  a linear  concentration  and 
density  gradient  and  a smooth  viscosity  gradient.  That  is,  the  profiles 
contained  some  steps. 

To  test  this,  the  results  of  runs  18  and  19  were  used.  In  these 
runs  a sharp  interface  existed  above  the  sink  elevation.  The  bottom 
layer  had  a viscosity  of  120  cp  while  the  top  layer  had  a viscosity  of 
4 cp.  If  buoyancy  forces  predominated,  flow  would  only  occur  below 
the  interface  (see  cloudy  layer  in  Figure  54).  Because  of  the  viscosity 
difference,  most  of  the  flow  occurred  above  the  interface. 

The  profile  resembles  those  obtained  from  the  one-dimensional, 
two-layer  case  in  Chapter  IV.  The  profiles  are  compared  in  Figure  55 
and  are  shown  to  be  very  close.  Note  that  the  profile  in  the  two-layer 
case  shows  a much  higher  peak  than  in  cases  with  a smooth  viscosity 
gradient.  It  was  concluded  that  the  peaks  in  the  observed  profiles 
were  due  to  steps  in  the  density  and  viscosity  profiles. 

So  while  the  correlations  between  the  observed  and  predicted 
velocity  profiles  are  reasonably  good,  the  results  would  be  even  better 
if  the  back  wall  were  farther  upstream  (or  removed  completely)  and  if 
steps  did  not  exist  in  some  gradients.  The  mathematical  model  proposed 
for  viscosity-stratified  flow  in  this  report  does  account  for  the 
important  trends  observed  in  the  experiments. 


125 


V 4 


Error  Analysis 


Three  variables,  B,  C,  and  u(y),  were  calculated  from  the  experi- 
mental data  and  used  in  verifying  the  mathematical  model.  The  results 
of  an  error  analysis,  performed  to  determine  the  effect  of  error  in 
experimental  measurements  on  the  model  results,  are  presented  below. 

The  error  in  the  variable  u can  be  expressed  as  6u.  If  u is  a 

function  of  variables  x, , x„,  ...  x with  error  6x, , 6x„,  6x  , then 

12  n 12  n 


<5u 


3u 

9x. 


6x. 


(110) 


The  error  contribution  of  the  ith  variable  is 


3u  . 

~ OX. 

3x  i 


The  percent  error  can  be  expressed  as 


6u 


x 100% 


The  error  analysis  was  first  carried  out  on  the  horizontal 
velocity,  u.  u was  determined  using  Ax/At  where  Ax  is  the  displacement 
in  time  At.  The  position  of  the  dye  streak  was  read  to  the  nearest 
tenth  of  a centimeter  and  the  clock  was  read  to  the  nearest  second. 

The  error  analysis  shown  in  Table  3 indicates  that  the  expected  maximum 
error  for  a typical  velocity,  assuming  that  roundoff  error  in  the 
parallax  correction  program  is  small,  is  9.1%.  Most  of  the  error 
originates  in  measuring  the  displacement  of  the  dye  streak. 

In  calculating  the  values  of  the  parameters  B and  C,  it  is 


necessary  to  measure  the  concentration  gradient  of  the  substance 
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TABLE  3 

ERROR  ANALYSIS  FOR  VELOCITY  (u) 


Parameter 

Typical 

Value 

Error 

% 

Error 

Error 

Contribution 

% Error 
Contribution 

Ax 

1.2 

0.1 

8.3 

0.00167 

8.3 

At 

60. 

0.5 

0.83 

0.00017 

0.8 

u 

0.02 

0.00184 

9.1 

causing  the  density  and  viscosity  gradient.  This  measurement  was  sub- 
ject to  error  in  the  instrument  used  to  determine  the  concentration  and 
in  exactly  measuring  the  elevation  at  which  the  sample  was  taken. 
Another  problem  is  that  the  gradient  was  treated  as  a constant;  yet  the 
fluid  was  placed  in  the  flume  in  discrete  layers.  These  layers  may  not 
have  faded  into  a smooth  gradient  by  the  time  of  the  run. 

The  confidence  limits  of  a slope  of  a regression  line,  where  b 
is  the  regression  coefficient  for  the  slope,  are 

b * Ca/2  S~  "Vn  - 1 
x ’ 

where 

tfl^2  “ student's  t for  significance  a/2 

- standard  error  of  estimate 
e 

■ standard  deviation  of  x 
n ■ number  of  observations 
For  this  study  a * 0.025. 

Using  the  above  formula  for  b ■ dy/dy  and  dp/dy,  the  error  in 
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the  viscosity  and  density  gradients  was  estimated.  The  error  in  de- 
termining dp/dy  was  virtually  always  less  than  30%  and  in  determining 
dp/dy  was  less  than  20%. 

The  accuracy  of  the  density  (specific  gravity)  measurement 
using  the  Pycnometer  method  according  to  ASTM  Standard  D 1429  (1971b) 
is  +0.005  specific  gravity.  The  viscometer  is  accurate  to  1%  while  the 
conductivity  meter  is  accurate  to  1%  according  to  the  manufacturer.  The 
less  accurate  of  the  ref ractometers  is  accurate  to  0.5%  but  it  was 
difficult  to  obtain  readings  with  better  than  1%  repeatability. 

Length  measurements  were  accurate  to  the  nearest  0.5  cm  while 
the  accuracy  of  the  flow  measure  was  taken  as  10%  since  it  was  determined 
by  integrating  the  velocity  profile.  Length  could  be  read  to  the  nearest 
0.1  cm  but  it  varied  during  the  run  and  was  recorded  to  the  nearest 
centimeter. 

The  results  of  the  error  analysis  for  the  parameter  C shown  in 
Table  4 indicate  that  most  of  the  error  is  contributed  by  the  measure- 
ment of  the  viscosity  gradient.  As  mentioned  earlier,  this  was  due  to 
the  fact  that  the  viscosity  gradient  did  not  always  fall  on  a straight 
line  on  a semi-log  graph. 

The  error  analysis  for  B is  given  in  Table  5.  The  principal 
source  of  error  was  the  inability  to  obtain  a perfectly  linear  density 
gradient.  The  characteristic  length  also  contributed  some  error  since 
it  varied  during  the  run  and  was  only  recorded  to  the  nearest  centimeter. 

While  the  percent  error  on  both  B and  C is  large,  it  must  be  re- 
membered that  these  values  represented  the  maximum  error,  and  almost  all 
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TABLE  4 


ERROR  ANALYSIS  FOR  C 


Parameter 

Typical 

Value 

Error 

% 

Error 

Error 

Contribution 

% Error 
Contribution 

l 

20. 

0.5 

2.5 

0.05 

2.5 

V 

2.00 

0.02 

1. 

0.02 

1. 

dp/  dy 

0.20 

0.06 

30. 

0.60 

30. 

C 

2.00 

0.67 

33. 

TABLE  5 

ERROR  ANALYSIS  FOR  B 

Parameter 

Typical 

Value 

Error 

% 

Error 

Error 

Contribution 

% Error 
Contribution 

l 

20. 

0.5 

2.5 

0.18E08 

10 

c 

0.0012 

0.00024 

20. 

0.38E08 

20 

Po 

1.00 

0.01 

1. 

0.18E07 

1 

po 

2.00 

0.02 

1. 

0.18E07 

1 

q 

0.05 

0.005 

10. 

0.18E08 

10 

B 

1.88E08 

0.78E08 

42 
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of  the  error  is  contributed  by  error  in  fitting  a slope  to  the  density 
and  viscosity  gradients.  In  general  the  error  was  much  less  than  the 

values  used  in  the  error  analysis.  i 

The  velocity  profiles  are  in  general  not  extremely  sensitive  to 
changes  in  B.  The  only  zone  in  which  B affects  the  profile  significantly 
is  the  transition  zone  from  buoyancy  to  viscosity  dominated  flow.  An 

/ 

error  in  B would  affect  the  x-coordinate  at  which  this  transition  occurs. 

The  velocity  profile  for  run  57  (Figure  53b)  appears  to  be  the  only 
profile  seriously  affected  by  this  type  of  error.  A 42%  maximum 
error  is  not  extremely  significant  when  one  realizes  that  in  the  experi- 
ment, B varied  over  eight  orders  of  magnitude.  The  effect  of  B on  skew- 
ness can  be  found  in  Figures  29  through  31  in  Chapter  VIII. 

Since  the  flume  was  filled  to  produce  a linear  concentration  and 
density  gradient,  a linear  or  log-linear  viscosity  gradient  could  not  be 
achieved  since  viscosity  does  not  vary  linearly  or  log-linearly  with 
concentration.  Therefore,  it  was  necessary  to  settle  for  whatever 
viscosity  gradient  resulted  from  the  linear  concentration  gradient. 

This  type  of  error  shows  up  most  dramatically  in  run  50.3  (Figure  53d). 

The  viscosity  gradients  in  most  runs  contained  much  less  error  than  the 
33%  shown  in  Table  4.  The  effect  of  C on  skewness  can  be  found  in 
Figures  13  and  14. 
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CHAPTER  XII 


EXTENSION  TO  WEIR  FLOW 


In  most  practical  problems  the  type  of  line  sink  that  is  of  most 
interest  is  a weir.  If  a stratified  fluid  is  flowing  over  a weir  with  a 
substance  whose  concentration  gradient  is  known  upstream,  the  problem  is 
then  to  predict  the  concentration  of  that  parameter  as  it  flows  over  the 
weir. 


Problem  Formulation 


This  problem  is  solved  using  the  log  viscosity  gradient  model 

K + c2  ifj  + fl  - 0 (111) 

3Y  3Y  3Y 


The  principal  differences  between  flow  to  a line  sink  at  mid-depth  and  a 
weir  are  the  boundary  conditions.  The  weir  is  merely  a line  sink  lo- 
cated at  the  surface  rather  than  at  mid-depth.  The  flow  is  also  allowed 
to  slip  along  the  top  boundary,  since  this  represents  an  air-water 
interface.  The  boundary  conditions,  shown  in  Figure  56,  can  be  written 
as 


no,o) 

= 0 

(112) 

nx,o) 

- -0.5  , 

X>0 

(113) 

no.Y) 

» +0.5  , 

Y<0 

(114) 

nx.D 

- +0.5  , 

X>0 

(115) 

A 
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ff(x,l)  = 0 

(116) 

^|(X,0)  = 0 

(117) 

The  first  condition  (equation  112)  represents  the  weir;  the  next  three 
(equations  113,  114,  and  115)  represent  no  flow  through  the  top,  bottom, 

or  downstream  boundary;  the  fifth  (equation  116)  represents  drag  on  the 
bottom;  and  the  sixth  (equation  117)  represents  slip  along  the  surface. 

Note  that  the  solution  presented  is  valid  only  for  the  case  in 
which  the  depth  of  flow  over  the  weir  is  small  in  comparison  to  the 
total  height  over  the  weir,  since  the  sink  in  the  model  is  treated  as 
being  only  one  grid  cell  wide.  A weir  with  large  head  can  be  modeled 
by  using  a sink  several  cells  wide,  provided  the  inertial  forces  are 
small. 


Velocity  Profiles 

The  equations  were  solved  numerically  using  the  same  techniques 

as  with  the  line  sink.  First  the  parameter  B was  varied  to  determine 

how  it  would  affect  the  flow.  The  streamlines  of  the  solution  for 
3 5 7 

B = 10  , 10  , 10  are  shown  in  Figure  57.  As  was  demonstrated  earlier, 
for  large  values  of  B,  the  withdrawal  layer  is  very  narrow,  since 
buoyancy  forces  predominate.  As  B decreases,  the  viscous  effects 
become  greater  and  the  profile  reaches  the  bottom  of  the  duct. 

Figure  58  gives  the  velocity  profiles  at  X ■ 10  for  the  stream- 
lines shown  in  Figure  57.  For  large  B the  velocity  tends  to  be  restrict- 
ed to  the  surface  layers. 

Next  the  parameter  C was  varied  (0.1,  2.0,  5.0)  for  a constant 


WEIR  MODEL 
C - 2.0  X - 10.0 


Fig.  58.  Effect  of  B on  velocity  profiles 
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value  of  B = 10^.  The  resulting  velocity  profiles  at  X = 10  are  given 
in  Figure  59.  The  effect  of  C on  these  profiles  is  not  significant 
until  C becomes  fairly  large  as  was  shown  in  the  one-dimensional  results 
in  Figure  10.  Since  X is  relatively  small  for  B = 10^ , buoyancy  forces 
tend  to  mask  the  effect  of  the  viscosity  stratification.  As  X becomes 
large  the  profiles  do  approach  those  given  in  Figure  10. 

Effluent  Concentration 

The  usual  problem  of  interest  in  selective  withdrawal  over  a 
weir  is  to  predict  the  concentration  of  some  substance  in  the  effluent 
knowing  the  concentration  gradient  upstream.  The  concentration  j.n  the 
effluent  can  be  given  by 

l 

s e = J u(y)s(y)dy  (118) 

o 

where 

se  = effluent  concentration 
s(y)  = concentration  at  cross  section  of  interest 
For  the  following  analysis,  dimensionless  concentration 
S = s/s (1/2)  will  vary  linearly  with  dimensionless  depth  as  given  by 
S = 2Y  (119) 

If  a uniform  velocity  profile  exists,  then  Sc  = 1.  If  the  velocity  is 
greater  near  the  surface,  then  the  effluent  concentration  will  be  less 
than  one.  Dimensionless  effluent  concentration  is  therefore 


U(Y)2YdY 


(120) 
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The  effluent  concentration  depends  on  the  distance  from  the  weir 
at  which  the  concentration  gradient  is  specified.  If  the  gradient  is 
specified  far  upstream  (X  -»■  <*>) , then  the  effluent  concentration  is  purely 
a function  of  C.  This  relationship  is  given  in  Figure  60,  using  the 
solution  from  Chapter  IV. 

If  the  gradient  is  specified  at  some  finite  value  of  X,  the 
effluent  concentration  depends  also  on  B.  This  case  represents  that  of 
a tank  with  a fixed  back  wall. 

The  effect  of  B on  the  effluent  concentration  is  given  in  Figure 
61.  When  B is  large  or  X is  small,  buoyancy  forces  tend  to  skim  layers 
of  low  concentration  off  the  top  of  the  fluid.  As  one  moves  farther 
from  the  weir,  the  effect  of  buoyancy  is  diminished.  For  X the 

effluent  concentration  is  only  a function  of  C.  As  was  shown  earlier, 
when  B is  small,  viscous  forces  dominate  at  smaller  values  of  X than 
when  B is  large. 


CHAPTER  XIII 


EXTENSION  TO  NON-NEWTONIAN  FLUIDS 


The  problem  of  stratified  flow  in  a non-Newtonian  fluid  is  more 
complicated  than  the  Newtonian  case  since  more  than  one  parameter  is  re- 
quired to  describe  the  relationship  between  shear  stress  and  shear  rate. 
When  the  flow  is  stratified,  all  of  these  parameters  may  vary  with  depth. 
In  this  chapter,  only  the  parameter  that  changes  dramatically  with  depth 
will  be  treated  as  a variable  in  modeling  the  flow.  A stratified  Bing- 
ham plastic  fluid  is  studied  for  two  cases:  yield  stress  varying  as  a 
function  of  the  logarithm  of  depth  and  yield  stress  varying  as  a 
linear  function  of  depth. 

For  a Bingham  plastic  fluid,  represented  by  the  model 


T = y -T—  + T , T >T 

yx  dy  oyx  ' 1 oyx 


(121) 


both  y and  can  vary  with  depth.  The  Bingham  plastic  model  is  impor- 
tant since  suspensions  of  fine-grained  sediments  behave  as  Bingham 
plastics. 

Concentration  of  fine-grained  suspended  solids  in  a dredged 

material  containment  area  have  been  reported  to  increase  with  depth  by 

Montgomery  (1978)  and  Walski  and  Schroeder  (1978) . This  would  imply 

that  both  y and  t Increase  with  depth. 

oyx  r 

In  Appendix  D,  it  is  shown  that  the  increased  suspended  solids 
with  depth  affect  the  yield  stress  more  than  the  viscosity  so  that  the 
y term  can  be  treated  as  a constant  in  the  model  developed  in  this 
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chapter.  This  means  that  once  the  yield  stress  Is  exceeded,  the  fluid 
behaves  as  an  unstratified  fluid. 

In  this  part  of  the  report,  a one-dimensional  model  is  developed 
for  a Bingham  plastic  fluid  in  which  yield  stress  varies  with  depth. 

The  model  predicts  depth  of  flowing  zone  and  velocity  profile  for  this 
type  of  flow. 


Logarithmic  Yield  Stress  Gradient 

The  definition  sketch  of  the  system  to  ’ e investigated  in  this 
chapter  is  given  in  Figure  62.  Note  that  the  boundary  conditions  have 
more  in  common  with  the  weir  model  than  the  one-dimensional  line  sink 
models.  The  yield  stress  varies  with  depth  as  given  by 


t = x y 
lyx 


oyx  ay 


The  model  to  be  investigated  can  be  given  by 

3x 


= !R 
3y  3x 


du 

T - — p -T-  + T 
yx  dy  oyx 


^ - 0 
dy 


with  boundary  conditions 

u(l)  = 0 

x (0)  = 0 
yx 


T > X 
yx  — oyx 


x < x 
yx  oyx 


I 


udy  » q 


These  equations  can  be  non-dlmenslonallzed  using  q,  1,  p as 


(122) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 


* •* 
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characteristic  flow,  length  and  viscosity.  The  non-dimensionalization 
gives  the  following  substitutions 


T,  ^ 

U = — u 

q 


*2r 


_32X 


qp 


*2X 


I = 


qw 


*2T 


oyx 


qp 


P = !_  iE 

x qp  3x 


9t  _ I?  3Tyx 
3Y  qp  3y 


Note  t without  subscript  yx  indicates  dimensionless  t 


yx 


x is  known  as 
a 


the  Bingham  number  at  Y = 1 and  is  a ratio  of  viscous  to  yield  force. 
Substituting 
3t 


3Y  ' Px 


(129) 


dU  , , , 

T = -rrr  + T , IT  > T 

dY  O 1 1 — O 


(130) 


dU  _ 
— = 0 
dY 


Boundary  conditions 

t (0)  =0 

U(l)  =■  0 
1 


/ 


UdY 


, | x | < T 


(131) 

(132) 

(133) 

(134) 


Treating  y as  a dimensionless  gives 

t (Y)  = t Yb 
o a 

Equation  129  can  be  solved  with  t(0)  • 0 to  give 

t = -P  Y 
x 


(135) 


(136) 
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There  are  two  cases,  the  first  is  where 

|t(Y)|  > tq  (137) 

where  the  fluid  behaves  as  a Newtonian  fluid  and 

|t(Y)|  < to  (138) 

where  the  velocicy  gradient  is  zero.  These  cases  give  rise  to  three 
types  of  flow 

1.  | x (Y) | tq  for  all  Y,  the  fluid  behaves  as  a Newtonian 
fluid  and  the  solution  given  in  Chapter  IV  is  valid. 

2.  | t (Y) | < Tq  for  all  Y,  the  fluid  behaves  as  a solid  plug. 

The  boundary  condition  U(l)  * 0 forces  the  plug  to  be  stationary. 

3.  There  exists  Y = L such  that 

| x (Y) | > tq  for  Y<L 

|t(Y) | < t for  Y>L 
o 

which  combined  with  the  boundary  condition  (U(l)  = 0),  means  that  there 
is  a stationary  zone  on  the  bottom  with  flow  above.  This  type  of  flow 
is  examined  below  in  greater  detail.  L is  referred  to  as  the  depth  of 
the  flowing  zone. 

L can  be  determined  by  noting  that  at  L 

|t(Y)|  = to  (139) 

which,  upon  substituting  from  equations  135  and  136,  implies 
* h 1 

' ~ ^ L , 0<L  < 1 (140) 

a 

For  Y>L 


37  - 0 (131) 

which  implies 

U(Y)  - 0 , Y>L  (141) 
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For  Y^L  the  Bingham  plastic  model  gives 


- § + x Yb  = P Y 
dY  a x 


which  can  be  solved  to  yield 


x L 
a 


b+1 


b + 1 


U = 


r - ■ 


P L 
x 


a - * 


L<Y<1 


(142) 


0<Y<L 


(143) 


The  continuity  boundary  condition  (equation  134)  can  be  used  to  solve 
for  P 

'C 

r 

, b+2 


- P = ~ (l  + T ~ 

X L3  ^ a b + 2 

. b-1 


(144) 


Recalling  - P^/x^  = L“  x,  the  above  relation  can  be  substituted  for  P^ 


and  solved  for  L 

L 


3_  fb+2\ 

Ta  H 


1/ (b+2) 


(145) 


The  above  equation  gives  the  depth  of  the  forward  flowing  layer  as  a 

function  of  x and  b only.  The  results  are  plotted  in  Figure  63. 

The  depth  of  the  flowing  2one  (L)  is  dependent  on  both  the 

yield  stress  of  the  fluid  on  the  bottom  (r  ) and  the  curvature  of  the 

relationship  between  yield  stress  and  depth  (b) . As  b approaches  one, 

it  becomes  less  likely  that  the  case  in  which  part  of  the  fluid  moves 

and  part  remains  stationary  will  occur.  This  type  of  zone  will  occur 

when  the  x and  x lines  Intersect  as  shown  in  Figure  62.  This  is 
yx  oyx 

much  more  likely  to  occur  when  the  is  highly  curved  (large  b) 

since  the  x^  line  is  straight. 

Velocity  profiles  for  the  case  in  which  xq  ■ are  shown  in 
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A- 
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t 

Figure  64  for  a series  of  values  of  b and  x . 

a 

Linear  Yield Stress  Gradient 


Another  potentially  important  case  is  that  in  which  there  exists 
a layer  of  water  (TQyX  = 0)  over  a material  with  a yield  stress  gradient 
which  increases  linearly  with  depth.  This  case,  which  represents  the 
type  of  flow  that  may  exist  in  a dredged  material  containment  area  with 
a sharp  interface  between  clear  water  and  a mud  layer,  is  shown  in 
Figure  65.  Given  that  the  interface  is  at  depth  L l,  the  yield  stress 

cl 

can  be  expressed  as 


oyx 


dT  (y-JlL  ) 
oyx  a 

dy 


y<£L 
— a 


Non-dimens ionalizing  as  earlier  yields 

0 , Y<L 


x = 
o 


a(Y-L  ) 


Y>L 
— a 


where 


a = dx  /dY 
o 


The  equations  to  be  solved,  in  dimensionless  form,  are 

- — * p 
3Y  x 


dU  . 

T “ dY  + To 


(L  <Y<L) 
a 


dU  „ 

= 0 

dY 


(Y>L) 


dU 

dY 


Y±La 


(146) 


(147) 


(129) 

(148) 

(149) 


(150) 


= 2.0 


Fig.  64.  Effect  of  Bingham  number  on  velocity  profiles 
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To  = a(Y_La)  » !lY>La  (151) 

with  boundary  conditions  given  by  equations  132,  133,  and  134.  t and  U 

are  continuous  at  L . 

a 

As  in  the  log  type  yield  stress  gradient,  three  types  of  flow 
may  exist 

1-  (L  = L)  Newtonian  flow  above  Y = L;  no  flow  below. 

a 

2.  (L  =*  1)  Newtonian  flow  throughout  fluid. 

3.  (L  >L>1)  forward  flowing  layer  penetrates  into  the  portion 
of  the  fluid  that  is  a Bingham  plastic. 

The  x-momentum  equation  can  be  solved  to  yield  (with  t(0)  = 0) 

t = -PxY  (136) 

L can  be  found  by  substituting  for  t in  the  above  equation  and  noting 

that  the  case  that  is  of  interest  is  L <L 

a— 

a(Y-L  ) = -P  Y for  Y = L (152) 

a X 

substituting  L for  Y yields 
aL  L 

L = a + P = 1 + (P  /a)  (153) 

x x 

For  Y>L 


dU 
dY  = 


0 


(149) 


which  implies  U(Y)  * 0 

For  Y<L,  there  are  two  possible  velocity  profiles  Y>L  , Y<L  . For 

H cl 

Y<L  , x * 0 and  r(L  ) ■ -P  L 
— a o 'a  x a 

The  last  condition  serves  as  a top  boundary  condition  for  the  case  Y>L 

Substituting  for  t using  the  Bingham  model  and  t = a(Y-L  ) 

o a 
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^ + t = - ^ + a(Y-L  ) = - P Y 
dY  o dY  a x 


(154) 


which  can  be  solved  using  U(L)  = 0 to  giye 


U - <px  + a)  T 


L <Y<L 
a 


i 


- a L L 
a 


i - 


Substituting  L into  equation  155  gives 
a 


” - <p*  + *)  f- 


L.) 2 


- a L L 7^-1 
a L 


(155) 


(156) 


which  is  the  boundary  condition  needed  to  solve  for  the  velocity  in  the 


Newtonian  depths  (Y<L  ).  That  solution  is 

a 


P L 
x a 


7-  I - 1|  + (Px  + a)  f- 


- 1 


- a L L 
a 


f L 

C 

lL 


- 1 


0<Y<L 
— a 


In  summary,  the  velocity  profile  can  be  written 


/ 


U = 


Ux  = 0 


U2  - <PX  + a)  f 


, Y>L 


P L 

U = 

3 2 


2 r 


Y_ 

L 


Y 

L 


- 1 


- a L L 
a 


(157) 


(158) 


+ U,(L) 
2 a 


l - • L«iYi-L 


0<Y<L 
— a 


To  eliminate  P^  from  the  equation,  the  continuity  boundary  condition  is 
used 


L L 

J U2dY  + J U. 


dY 


(159) 
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The  integrals  can  be  evaluated  to  give 


(160) 


The  expression  for  can  now  be  substituted  into  equation  153  and 
solved  for  L 


L • VLa  + IT  <161> 

a 

The  above  relationship  can  be  plotted  as  shown  in  Figure  66  which  can 

be  used  to  determine  the  depth  of  the  flowing  zone,  knowing  only  a and 

L . 
a 

Velocity  profiles  generated  using  the  above  model  for  L =0.2 

a 

are  shown  in  Figure  67 . As  a becomes  larger  the  velocity  profiles  are 
cut  off  much  sharper  below  the  interface.  Figure  68  shows  that  the 
velocity  profiles  may  not  be  identical  even  though  the  withdrawal  depth 
and  flow  are  the  same. 


Implications 


The  above  solutions  can  be  used  to  describe  flow  in  a dredged 
material  containment  area  at  locations  far  from  the  weir  and  at  times 
when  wind-generated  currents  are  negligible.  Once  suspended  solids  in 
the  area  have  settled  to  depths  greater  than  L,  they  are  removed  from 
the  water  column.  Unless  they  are  resuspended  by  wind-generated  tur- 
bulence, they  will  not  flow  over  the  weir. 

The  velocity  profile  can  be  predicted  knowing  only  the  flow, 
dimensions  of  the  basin,  concentration  gradient,  and  relationship  of 
concentration  to  yield  stress.  The  log  gradient  case  should  be 


V * 


Fig.  66.  Depth  of  flowing  layer  vs.  gradient 
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V * 


Fig.  67.  Velocity  profiles  for  L = 0.2 
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applicable  to  areas  without  a sharp  Interface  while  the  linear  gradient 
should  be  appropriate  when  a sharp  interface  between  a clear  water 
layer  and  a layer  with  high  solids  concentrations  exists.  If  the 
yield  stress  gradient  does  not  fit  one  of  the  cases  solved  in  this 
chapter  of  the  report,  the  equations  of  motion  can  be  solved  with  the 
appropriate  yield  stress  gradient. 

In  the  area  near  the  weir,  buoyancy  forces  predominate  and  the 
model  presented  in  the  previous  chapter  should  be  appropriate.  The 
model  presented  in  this  chapter  is  appropriate  when  buoyancy  and 
inertial  forces  are  no  longer  important.  This  will  occur  when  the 
streamlines  are  horizontal.  The  location  at  which  this  occurs  is 
given  in  Chapter  XII. 


CHAPTER  XIV 


CONCLUSIONS  AND  OBSERVATIONS 


Velocity  profiles  in  a fluid  in  which  the  viscosity  varies 
vertically  will  be  different  than  those  for  a fluid  with  constant 
viscosity.  This  vertical  variation  in  viscosity  is  referred  to  as 
viscosity  stratification. 

In  regions  of  a fluid  where  buoyancy  and  inertial  forces  are 
of  a smaller  magnitude  than  viscous  forces,  viscosity  stratification 
will  cause  the  velocity  profile  to  skew  positively  in  the  direction 
of  low  viscosity. 

When  flow  is  to  a line  sink  located  at  mid-depth  in  a duct,  the 
skewness  can  be  represented  by  the  third  moment  of  the  dimensionless 
velocity  distribution  about  the  sink  elevation.  This  unarbitrary 
indicator  of  skewness  can  be  defined  as 


skew 


L/2 


U(Y)Y  dY 


-1/2 


r 


U(Y)dY 


-1/2 


(57) 


For  viscosity-stratified  flows  with  Reynolds  number  less  than 
100,  the  flow  can  be  modeled  by  the  dimensionless  equation 


1L.  J JL  ill  + fll  *2v\ + 8*  epo  il 

Wo  3Y4  \.Wo  dyJ  3Y3  9y 2 J 3Y2  ^o  8X 


(26) 
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where  the  first  term  represents  viscous  forces,  the  next  two  represent 

viscosity  stratification  and  the  final  term  represents  buoyancy  forces. 

For  significantly  large  values  of  the  A/pq  3y/3y  parameter, 

viscosity  stratification  will  occur  at  some  distance  from  the  sink. 

The  distance  at  which  viscosity  stratification  becomes  significant  is 

4 

dependent  on  the  parameter  B • g l zp^/qv^.  That  is,  buoyancy  forces 
may  mask  the  viscosity  stratification  near  the  sink. 

In  viscosity-stratified  flow  there  exists  a region  in  which 
viscous  forces  dominate  at  lower  elevations  in  the  fluid  while  buoyancy 
forces  dominate  in  the  top  layers.  This  type  of  flow  is  referred  to 
as  having  "incipient  skewness." 

The  velocity  profiles  in  a two-dimensional  duct  approach  the 
profiles  for  one-dimensional  flow  as  the  distance  from  the  sink  increases. 
Skewness  can  develop  in  a flow  even  with  a fixed  back  wall. 

The  above  conclusions  were  made  using  results  from  both  the 
mathematical  model  and  the  experimental  model. 

With  the  appropriate  boundary  conditions,  the  model  for 
viscosity-stratified  flow  can  be  used  to  predict  velocity  profiles  for 
flow  over  weirs,  when  the  conditions  assumed  in  developing  the  model 
exist. 

Velocity  profiles  in  a Bingham  plastic  with  varying  yield 
stress  can  be  dramatically  affected  by  this  variation  in  yield  stress. 

A mathematical  model  for  this  type  of  flow  may  be  used  to  predict 
motion  of  cohesive  sediments.  Algebraic  equations  were  developed 
relating  yield  stress  gradient  to  the  depth  of  the  flowing  zone. 
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PROGRAMS  FOR  NUMERICAL  SOLUTIONS 


APPENDIX 


PROGRAMS  FOR  NUMERICAL  SOLUTIONS 


The  equations  presented  earlier  in  the  text  were  solved  using  a 
series  of  computer  programs  each  with  a different  combination  of 
boundary  conditions,  viscosity  gradients  and  simplifying  assumptions. 

The  programs  all  had  the  same  overall  layout  shown  in  the  flowchart  in 
Figures  69  and  70. 

The  program  first  reads  the  data  for  the  run,  echoes  the  data, 
and  calculates  initial  and  boundary  conditions  and  other  constants  for 
the  solution.  It  then  prompts  the  user  with  a question  mark.  The  user 
then  selects  whether  to  advance  the  solution  through  N iterations;  output 
stream  function,  dimensionless  velocity  profile,  scaled  velocity  profile, 
or  skewness  at  any  combination  of  cross  section;  output  the  coordinates 
of  the  stream  function;  or  stop  the  program. 

A list  of  the  programs  developed  and  their  features  is  given 

below. 

DUCT  - constant  density;  linear  upstream  velocity  profile  and 
viscosity  gradient;  inertia  terms  included. 

TANKU  - constant  density;  fixed  back  wall;  linear  viscosity 
gradient. 

BUOY  - variable  density  included;  linear  viscosity  gradient; 
slip  at  boundary. 

DRAG  - same  as  BUOY  with  no  slip  at  boundary. 

TRY  - same  as  DRAG  with  log  viscosity  gradient. 

VISC  - same  as  TRY  except  y = f(’f)  instead  of  f(Y). 
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Fig.  69.  Flowchart  for  numerical  solutions  - 1 


r + 
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Fig.  70. 


Flowchart  for  numerical  solutions  - 2 
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There  were  a 
dimensional  problems 


number  of  smaller  programs  written  to  solve  one- 
simplifying  the  evaluation  of  long  expressions. 


These  are  too  trivial  to  merit  listing  in  the  report 


APPENDIX  B:  VELOCITY  DATA  REDUCTION 


In  determining  the  location  of  the  dye  streak  from  the  photo- 
graph, a correction  must  be  made  due  to  the  fact  that  light  bends  in 
passing  from  one  medium  to  another  and  that  there  is  a finite  distance 
between  the  dye  streak  and  the  back  wall  of  the  flume  on  which  the  grid 
is  located.  In  this  appendix,  a formula  for  correcting  the  apparent 
position  of  the  dye  streak  is  developed. 

The  apparent  location  of  the  dye  streak  as  given  on  the  slide 
is  designated  (x^,  y ) where  x and  y are  horizontal  and  vertical  dis- 
tances, respectively,  from  the  point  perpendicular  to  the  camera  (the 
camera  location  is  fixed  for  each  flume  (Figure  71)).  The  distance  from 
this  center  point  dA  is  given  by 

dA  " V\  + yI  <162> 

Since  the  thickness  of  the  glass  or  plexiglass  is  very  small 
in  relation  to  the  thickness  of  water  and  air  through  which  the  light 
must  pass  and  the  index  of  refraction  of  glass  and  plexiglass  is  not 
greatly  different  from  water,  the  refractance  due  to  the  tank  walls 
will  be  neglected.  Furthermore,  the  sugar  and  CMC  have  a negligible 
effect  on  the  refractance  of  water  compared  to  the  differences  between 
the  refractance  of  air  and  water. 

Consider  Figure  72  in  which  x^  « 0 so  that  d^  ■ y^.  It  is 
desired  to  find  y,  the  actual  location  of  the  dye  streak  in  the  center 
of  the  tank  (w/2  from  back  wall). 
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Fig.  71.  Definition  sketch  (x-y  plane) 
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Snell's  Law  of  Refraction  states  that  the  ratio  of  the  sine  of 
the  angle  of  incidence  (sin  0 ) to  the  angle  of  refraction  (sin  By)  is 
inversely  proportional  to  the  absolute  indices  of  refraction  of  A (air) 
and  W (water) . 

Sln  9 A = 1.333  (163) 

sin  9„  1.000 

w 

6 and  0..  are  shown  in  Figure  72. 

A W 

0TT  is  therefore 
W 

0..  = arc  sin 

w 

Rules  of  similar  triangles  can  be  used  to  form  the  equation 


7 = yA  " 1,2  W t3n  9W 

(164) 

yA  “ y 

from  tan  0W  - w/2 

(165) 

tan  9A  = L 

(166) 

When  W<<L,  it  is  possible  to  replace  0^  by  ©c  where 

ec  = arc 

(167) 

0tT  can  then  be  written 

W 

r sin  ©J) 

6W  = arc  sinl  1>333  J 

(168) 

given  a series  of  equations  from  which  y can  be  calculated 

given  yA,  L, 

and  W. 

Since  in  most  cases  x ^ x^,  xA  must  also  be  corrected.  Defining 
ip  = arc  tan  (y^/x^)  (169) 


sin  0, 


1.333 


(164) 


it  is  possible  to  solve  the  above  equations  for  d knowing  d^  instead  of 


x knowing  and  to  convert  the  results  using 
x = d cos  <f> 
y = d sin 


(170) 

(171) 


The  working  set  of  equations  to  obtain  (x,  y)  given  (x  , y ), 

A A 


L,  and  W is 
d 


-V* 


XA  + yA 


4>  = arc  tan  (yA/ x^) 

0C  = arc  tan  (d^/L+W) 

0W  = arc  sin  (sin  0C/1.333) 


d = dA  - 1/2W  tan  0W 


(162) 

(169) 

(167) 

(168) 
(165) 


x = d cos  (170) 

y = d sin  <f>  (171) 

Given  the  location  of  the  velocity  profile  at  any  two  times,  it 
is  a complicated  and  tedious  process  to  correct  from  parallax  error 
and  then  using  the  correction  locations  to  calculate  velocity  profiles. 
It  is  even  more  difficult  when  the  profile  is  made  by  averaging  several 
profiles  each  of  which  must  be  done  separately. 

To  circumvent  these  problems,  the  program  VELPRO  was  prepared. 
VELPRO  consists  of  a calling  program,  VELPRO,  which  reads  the  data  and 
four  subroutines : 

1.  PARA  which  performs  parallax  error  corrections  for  each  point 
input ; 

2.  INTER  which  interpolates  between  the  points  to  form  corrected 


dye  streaks; 
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3.  PRO  which  calculates  the  velocity  between  each  combination  of 
two  streaks,  and  averages  these  values  to  obtain  an  average  maximum  and 
minimum  velocity  at  each  point; 

4.  STAI  which  calculates  the  skew  of  the  average  profile  and 
integrates  the  profile  to  obtain  the  flow  rate. 


PROPERTIES  OF  MODEL  FLUIDS 


This  Appendix  contains  data  on  the  concentrative  properties  of 
the  sugar  and  CMC  used  in  the  experimental  model.  The  procedure  for 
determining  viscosity,  density,  refractance,  and  conductivity  are  pre- 
sented, followed  by  a check  on  the  Newtonian  behavior  of  7H  CMC. 

Viscosity  Procedure 

The  viscosity  was  measured  using  Brookfield  Synchro-Lectric 
viscometer  models  LVF  and  LVT,  as  pictured  in  Figure  73.  Tests  were 
run  using  spindle  Nos.  1 and  4,  and  the  UL  Adaptor. 

Spindle  4 was  used  for  high  viscosities;  spindle  1 was  used 
for  intermediate  viscosities;  and  the  UL  adaptor  was  used  for  \ is- 
cosities  less  than  10  cp.  Spindles  1 and  4 are  shown  in  the 
rack  in  the  front  right  in  Figure  73,  whereas  the  sleeve  for  the  UL  is 
shown  on  the  viscometer  in  the  figure. 

The  readings  were  taken  by  lowering  the  viscometer  into  the 
sample  of  known  concentration  to  the  designated  level.  Readings  were 
taken  at  as  many  rotational  speeds  as  produced  readings  away  from  the 
ends  of  the  scale,  as  repeatability  could  not  be  achieved  at  the  low 
end  of  the  scale. 

The  temperature  was  measured  and  the  relative  viscosity  was 
obtained  by  dividing  the  viscosity  by  the  viscosity  of  water  at  that 
temperature.  The  relative  viscosity  was  plotted  vs.  concentration. 

During  the  viscosity  tests,  a sample  of  National  Bureau  of 
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Standards  oil,  standard  number  S20,  was  tested.  Since  the  volume  of 
the  sample  (1  pint)  was  small,  it  could  only  be  tested  using  the  UL 
Adaptor.  The  oil  at  26°C  was  found  to  have  a viscosity  of  28.9  cp, 
which  agreed  well  with  the  given  value  of  29.2  cp. 

Density  Procedure 

The  specific  gravity  of  all  samples  was  read  using  an  ASTM 
Hydrometer  No.  151H-62  (ASTM  E 100-66,  1971a),  which  corresponds  to 
ANSI  standard  Z102. 1-1967.  The  results  for  sugar  agreed  well  with 
values  for  sucrose  as  reported  in  the  CRC  Handbook  of  Chemistry  and 
Physics  (Weast  1968) . 

Because  of  the  small  concentrations  of  CMC  used,  very  poor 
reproducibility  was  achieved  with  the  hydrometer  for  CMC.  The  tests 
were  repeated  using  a pycnometer  meeting  ASTM  D 1429-60  (1971b). 

These  results  were  used  in  the  study.  The  results  for  specific  gravity 
were  converted  to  density  when  the  temperature  of  the  samples  was  known. 

Refractance  Procedure 

Refractance  was  used  to  Indicate  the  concentration  of  sugar  in 
the  samples.  In  low  concentrations,  the  refractance  was  measured  using 
an  American  Optical  Goldberg  Temperature  Compensated  (T/C)  Ref ractometer , 
Model  10423.  The  readings  were  an  index  of  refraction  with  full  scale 
from  1.3330  to  1.3730. 

For  sugar  concentrations  greater  than  30%  by  weight,  an 
American  Optical  T/C  Hand  Ref ractometer , Model  10431  was  used.  This 
instrument  was  less  accurate  (±0.25°Brix)  but  had  a range  of  0-50 
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Brlx.  The  refractance  readings  agreed  well  with  values  reported  in  the 
CRC  Handbook  (Weast  1968). 

Conductivity 

The  concentration  of  CMC  solutions  could  be  determined  by 
measuring  the  conductivity  of  the  solutions  using  a Yellow  Springs 
Instrument  Co.  Model  31  Conductivity  Bridge.  The  probe  was  lowered  into 
the  sample,  the  instrument  dial  was  adjusted  to  maximize  the  "shadow" 
and  was  then  read  in  pmho/cm. 

No  temperature  correction  was  done  for  the  conductivity  since 
the  meter  and  flume  were  used  in  the  same  room  which  was  kept  at  25°C 
continuously. 

Results 

Figures  74,  75,  76,  77,  and  78  give  relationships  between  con- 
centration, specific  gravity,  relative  viscosity,  and  conductivity 
(at  25°C)  or  refractance  for  sugar,  7M,  7H,  7H4,  and  D-70  CMC.  In 
general,  the  refractance,  specific  gravity,  and  conductivity  were 
linear  functions  of  concentration  while  the  viscosity  increased  at  an 
increasing  rate  with  concentration. 

Non-Newtonian  Behavior 

An  investigation  was  made  of  the  7H  CMC  to  check  if  it 
exhibited  Newtonian  behavior  in  the  concentrations  used  in  this  study. 

It  was  felt  that  if  the  CMC  was  non-Newtonian,  it  would  obey  a power 
law  model. 


T — * 


Fig.  74 
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Fig.  76.  Properties  of  7H  CMC 
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The  rotational  speed  (ft)  in  radians/sec  and  shear  stress  at 
2 

the  bob  (t^)  in  degrees/cm  were  measured  for  an  array  of  rotational 
speeds,  cup  sizes  and  bob  sizes.  Note  that  the  bob  is  the  part  of 
the  viscometer  that  rotates  while  the  cup  is  the  container  holding 
the  fluid.  The  data  were  reduced  using  the  following  method  adopted 
from  Van  Wazer  et  al.  (1963). 

In  a rotational  viscometer  the  power  law  model  can  be  given  as 
xr0  = m|-du/dr|n  (172) 

The  shear  stress  can  be  related  to  the  torque  of  the  viscometer  by 


where  T is  torque,  dyne-cm 

h is  effective  length  of  spindle  with  end  corrections 
and  shear  rate  can  be  given  as 
du  dw 

= r 

dr  dr 

where  w is  angular  velocity  in  radian/sec. 

Substituting  for  du/dr  in  equation  174,  and  rearranging  gives 


(173) 


(174) 


, 1/n 


-dw 


2ir  hm 


r-(1+2">  dr 


(175) 


Integrating  with  w = 0 at  r = (cup  reading)  and  w e ft  at  (bob 
radius)  gives 

f t V/n/ 


2tt  hm 


J 


n \ L-2/n  p-2/n 

h “ c 


') 


(176) 


Taking  log  of  both  sides  and  substituting 


T - rb2irRgh  for  T 


(177) 
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where  x,  is  shear  at  the  bob 
b 

lnfi  = 1/n  ln(xb)  + In  1 - (RB/RC)2/n  (178) 

Plotting  lnfi  vs.  lnx^  allows  n to  be  determined  from  the  slope  of  the 
plots.  A 5 g/J l solution  of  7H  type  CMC  was  tested  for  a variety  of 
rotational  speeds  and  bob  sizes.  The  data  are  shown  in  Figure  79. 

A curve  fit  program  was  used  to  determine  the  slope  and  y- 
intercept  for  each  of  the  lines,  n was  found  to  be  1.01,  1.04  and 
1.01.  Since  these  values  of  n were  very  close  to  1,  n will  be  consid- 
ered 1 for  7H  CMC. 

m can  be  found  from  the  y-intercept  (y^) 

— y 

m = f ° [1  - (R  /R  )2]  for  n = 1 (179) 

z a c 

Since  n = 1,  m = p * actual  viscosity  v ■ 2.69  @28°C 

The  above  value  agrees  with  the  values  which  would  be  obtained  if  the 

fluid  were  assumed  to  be  Newtonian.  When  the  value  n is  one,  the  power 

law  model  reduces  to  Newton’s  Law  of  Viscosity.  It  was  concluded  that 

7H  CMC  was  Newtonian  at  5 g/£  or  lower  concentrations  for  shear 

2 

stresses  less  than  625  dynes/cm  (full  scale  on  viscometer). 

A plot  similar  to  Figure  79  was  prepared  for  7H4  CMC.  It  also 
proved  to  be  Newtonian  at  5 g/i. 


APPENDIX  D:  EFFECT  OF  SUSPENDED  SOLIDS  ON  VISCOSITY 


The  viscosity  of  a suspension  increases  with  an  increase  in  the 
suspended  solids  concentration.  This  increase  is  a consequence  of  the 
perturbing  of  the  flow  of  the  liquid  phase  around  the  solid  particle. 
The  suspended  particles  can  be  viewed  as  inhibiting  the  transfer  of 
momentum  through  the  fluid. 

Frisch  and  Simha  (1956),  Sardon  (1953)  and  Jeffrey  and  Acrivos 
(1976)  present  numerous  models  to  describe  the  dependence  of  viscosity 
on  concentration.  Th_  models  generally  treat  the  liquid  phase  as  a 
continuum  as  they  are  based  on  the  assumption  of  high  values  of  the 
Knudsen  number  (d/X,  where  d is  the  linear  dimension  of  the  apparatus 
and  X is  the  mean  free  path) . 

The  models  can  be  divided  into  three  categories  based  on  the 
assumptions  made  in  their  derivation.  The  first  category  is  "extremely 
dilute  suspensions"  in  which  the  suspended  particles  do  r.ot  interact 
with  one  another;  the  second  is  "dilute  suspensions"  in  which  there 
is  some  hydrodynamic  interaction  between  particles;  and  the  third  is 
"concentrated  suspension"  in  which  reactions  between  particles  become 
of  primary  importance. 

In  extremely  dilute  suspensions,  the  change  in  viscosity  due 


to  suspended  particles  is  the  sum  of  the  changes  caused  by  each  particle. 
The  concentration  dependence  of  viscosity  is  linear  in  this  range  and 
can  be  given  by  the  "intrinsic  viscosity,"  [p]  where 
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r i 1 dy 

[y  ] = — -j-  c -*■  o 
y ds 
o 


(180) 


[ y ] = intrinsic  viscosity,  cm  /gr 

y = viscosity  of  fluid  phase  without  suspended  particles, 

° 2 
dynes-sec/cm 

3 

s <=  concentration,  gr/cm 

The  models  of  extremely  dilute  suspensions  can  be  developed  by  linear- 
izing the  Navier-Stokes  equation  using  Stokes'  or  Oseen's  approximation. 
Since  there  is  no  interaction  between  particles,  the  boundary  condition 
on  flow  at  an  infinite  distance  from  the  particles  can  be  specified. 
Assuming  the  particles  to  be  non-interacting  spheres,  one  can  arrive 
at  Einstein's  solution 


Prel  = {r  = 1 + 2-5*  (181) 

o 

where  <j>  » volume  concentration  of  particles. 

The  solution  is  somewhat  different  if  one  assumes  dumbbell-shaped  par- 
ticles with  length,  L,  and  sphere  radius,  a. 

yrel  = 1 + 1.5(L/a)2<j>  (182) 

A wide  variety  of  relationships  have  been  derived  based  on  the  assump- 
tions used  in  constructing  the  models. 

When  the  solid  concentration  becomes  larger,  the  particles  begin 
to  exert  hydrodynamic  forces  on  one  another.  The  particles  will  also 
begin  to  aggregate.  This  can  be  seen  by  the  plotting  relative  vis- 
cosity vs.  concentration  as  done  by  Roscoe  (1953).  At  some  point,  the 
slope  begins  to  increase  with  concentration.  Models  to  describe  the 
behavior  of  particles  in  this  range  have  the  form 
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I 

I 

1 


i 


l'rel  - 1 + 2.5*  + X*2  (183) 

where 

14.1  Gold 
12.6  Simha 
7 . 35  Vand 
7.17  Vand 
4.51  Vand 
3.1  Vand 
2.5  de  Bruijn 

depending  on  the  researcher.  Where  aggregation  becomes  important,  the 
concentration  dependence  of  viscosity  can  be  written  as 

- 1 + s tv]  (1  + Kx  [pi  s + k2  Ip)2  s2  + ...)  (184) 

3 

where  s = weight  concentration,  g/cm  . 

Researchers  have  found  to  be  on  the  order  of  0.33  to  1.25.  The  last 
equation  can  be  rearranged  in  terms  of  "relative  specific  viscosity," 
p/s 

■ [p]  (1  + tv)  s + K2  [p ] 2 s2  + ...)  (185) 

where  p = specific  viscosity  - (p-pQ)/po 

Once  the  suspensions  become  highly  concentrated,  rational  models 
can  no  longer  be  developed  to  describe  the  viscosity  because  of  the 
complexity  of  interparticle  Interactions.  The  simplest  models  are  of 
the  form 


(186) 
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Since  suspensions  tend  to  become  non-Newtonian  in  this  range,  this  type 
of  model  is  of  very  little  value. 

Empirical  studies  by  Krone  (1963a,  1963b)  and  Migniot  (1968) 
have  shown  that  naturally  occurring  suspensions  exhibit  the  behavior 
of  Bingham  plastics  at  high  concentrations.  At  extremely  high  concen- 
trations, (>500  g/£)»  they  may  even  become  thixotrophic. 

Several  of  the  models  relating  concentrations  and  viscosity  are 
compared  with  field  data  gathered  by  Migniot  (1968),  Tiederman  and  Reisch- 
man  (1973),  and  Krone  (1963b).  These  values  are  shown  in  Table  6. 

Mignoit's  data  are  for  a variety  of  sediments  in  estuaries  and  rivers. 

The  data  of  Tiederman  and  Reischman  are  for  dredged  material.  Krone's  data 
are  for  San  Francisco  Bay  sediments.  These  experimental  values  are  compared 
with  values  calculated  using  models  for  extremely  dilute,  dilute,  and  con- 
centrated suspensions.  The  models,  taken  from  Frisch  and  Simha  (1956),  were 
extremely  dilute 


Vi  =p  --1*1  + 2. 5$ 
sp  rel  Y 

dilute 

W = 2.5<j>  + 14<|>2 
sp 

concentrated 
P®p  (l-4>)2,5 


(181) 


(187) 


(188) 


It  was  assumed  that  the  specific  gravity  of  the  particles  was  2.67. 

The  theoretical  models  consistently  predict  extremely  low 
specific  viscosities.  The  models  generally  Ignore  electrical  attraction 
and  assume  regular  shape  particles  which  is  generally  not  the  case  in 
naturally  occurring  sediments.  The  low  values  reported  by  Migniot  were 
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for  powdered  limestone.  There  is  fairly  good  agreement  between  pre- 
dicted viscosities  and  Mlgniot's  low  values. 

As  the  concentration  of  the  suspension  increases,  the  viscosity 
becomes  non-Newtonian.  This  will  occur  at  different  concentrations  for 
different  suspensions.  The  non-Newtonian  model  that  best  fits  the  sus- 
pension is  the  Bingham  plastic  model 


T * 

du 

— U — + T , 

|t  |>t 

(124) 

yx 

dy  y 

1 yx1  oyx 

du  * 

o 

It | <x  ^ 

(125) 

dy 

1 yx ' oyx 

where 


t » shear  force  in  x direction  in  plane  perpendicular  to  y, 
yX  dyne/ cm ^ 

u * velocity  in  x direction,  cm/sec 

2 

p * viscosity  coefficient,  dyne-sec/cm 

2 

Toyx  **  shear  stress,  dyne/cm 

The  yield  shear  stress  (which  Migniot  calls  "initial  rigidity") 
is  a function  of  material  type  and  concentration  and  is  an  indicator  of 
the  plasticity  of  the  suspension.  Migniot  (1968)  gives  the  yield 
shear  as 

x * n (189) 

o 

where 


s « concentration,  g/l 
n « factor  depending  on  material  type 
According  to  this  model  every  suspension  has  an  initial  rigidity  but 
many  are  very  low.  The  relationship  is  graphed  in  Figure  80,  taken 
from  Migniot  (1968).  Note  that  the  yield  stress  is  highest  for 
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Fig.  80.  Yield  stress  vs.  suspended  solids  concentration 
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kaolinite  and  lowest  for  powdered  limestone. 

Fine-grained  sediments  can  affect  the  behavior  of  Bingham 

plastics  in  both  the  y du/dy  term  and  the  toyx  term  in  the  Bingham 

model.  If  concentration  varies  linearly  with  depth,  it  is  clear  that 

the  t term  is  affected  much  more  than  the  y du/dy  term  since  y would 
oyx 

increase  linearly  with  depth  while  t would  increase  as  the  fourth 

oyx 

power  of  depth.  Therefore,  the  variation  of  viscosity  with  depth  in 
the  model  can  be  ignored  in  cases  where  the  Bingham  plastic  model 
applies.  This  is  done  in  the  models  presented  in  Chapter  XIII. 


APPENDIX  E:  NOTATION 


constants  in  one-dimensional  solution  with  log  viscosity 
gradient,  0 

— -jp-  , dimensionless  linear  viscosity  gradient,  0 

% ^ 

regression  coefficient  for  slope 

.4 

8*  ep0 

, dimensionless  ratio  of  buoyancy  to  viscous 

^ Wo  forces,  0 

slope  of  log  viscosity  gradient,  1/L 

c^  1,  dimensionless  log  viscosity  gradient,  0 

actual  distance  from  point  to  position  perpendicular  to 
camera,  cm 

apparent  value  of  d,  cm 

2 

diffusion  coefficient,  L /T 

constant  in  solution  by  Gehar  and  Mascalo,  1/LT 

value  of  function  f in  jth  interval  in  skewness  calcula- 
tion 

q/i‘  , dens ime trie  Froude  number,  0 

2 2 
980  cm/sec  , acceleration  due  to  gravity,  cm/sec 

length  of  viscometer  bob,  cm 

index  on  grid  cell  in  x-dlrection,  0 

index  on  grid  cell  in  y-dlrectlon,  0 

characteristic  length,  cm 

distance  from  camera  to  flume,  cm 

dimensionless  depth  of  forward  flowing  layer  (non- 
Newtonian)  , 0 
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dimensionless  depth  to  interface,  0 

number  of  grid  cells  in  x-direction,  0 

constant  in  power  law  model,  M/LT 

number  of  grid  cells  in  y-direction,  0 

exponent  in  power  law  model,  0 
2 

pressure,  M/LT 

2 2 

3p/3x,  pressure  gradient,  M/L  T 
3 

Pxi  /qyQ,  dimensionless  pressure  gradient,  0 

2 

characteristic  flow  rate,  cm  /sec 

2 

flow  at  cross-sections  x in  tank,  cm  /sec 

radial  distance,  cm 

qpo/viQ,  Reynolds  number,  0 

radius  of  viscometer  bob,  cm 

radius  of  viscometer  cup,  cm 
3 

concentration,  M/L 

3 

effluent  concentration,  M/L 

dimensionless  concentration,  s(Y)/s(l/2),  0 

dimensionless  effluent  concentration,  0 

standard  error  of  the  estimate,  0 

standard  deviation,  0 

time  to  reach  steady-state,  sec 

viscometer  torque,  dyne-cm 

velocity  in  x-direction,  cm/sec 

ui/q,  dimensionless  velocity  in  x-direction,  0 

dimensionless  velocity  in  jth  grid  cell,  0 

dimensionless  velocity  far  from  sink,  0 
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v 

V 

V* 

w 

U 


x 


X 

V 


velocity  in  y-dlrection,  cm/sec 

vi/q,  dimensionless  velocity  in  y-direction,  0 

dimensionless  velocity  at  which  tank  surface  drops,  0 

angular  velocity,  radians /sec 

width  of  flume,  cm 

horizontal  distance,  cm 

apparent  x location  of  point,  cm 

distance  to  back  wall  of  tank,  cm 

x/i,  dimensionless  horizontal  distance,  0 

x/x  , 0 
max 

constant  in  formula  for  viscosity  of  suspension,  0 

dimensionless  distance  to  back  wall  of  tank,  0 

vertical  distance,  cm 

Intercept  of  ft  vs.  plot,  0 

apparent  location  of  point,  cm 

y/l,  dimensionless  vertical  distance,  0 

point  about  which  skew  is  calculated,  0 

dimensionless  elevation  of  jth  grid  cell,  0 

thickness  of  withdrawal  zone,  cm 

l/pQ  dpg/dy,  density  gradient,  1/cm 

angle  of  Incidence  from  air,  radians 

corrected  angle  of  incidence  from  air,  radians 

angle  of  refractance  in  water,  radians 

mean  free  path  of  molecules,  L 

viscosity,  M/LT 


characteristic  viscosity  or  viscosity  of  fluid  phase  with 
no  suspended  particles  (Appendix  D) , M/LT 

viscosity  of  bottom  layer,  M/LT 

Ap/pq,  0 

viscosity  of  top  layer,  M/LT 
(p-po)/vo>  specific  viscosity,  0 
p/pQ,  relative  viscosity,  0 

3 

1/u  dp/dc  intrinsic  viscosity,  cm  /gr 
o 

3 

density,  gr/cm 

3 

characteristic  density,  gr/cm 

3 

density  of  fluid  with  no  motion,  gr/cm 

2 

l t /qp  , dimensionless  shear  stress,  0 
yx 

2 

i.  T0yX/<IW  » dimensionless  yield  stress,  0 

2 

shear  stress  at  viscometer  bob,  dyne/ cm 

component  of  shear  stress  tensor  (also  viscous  transport 
of  x-momentum  in  y direction),  M/LT^ 

2 

yield  stress,  M/LT 

arc  tan  (y/x),  parallax  error  correction,  radius 
volume  fraction  of  suspended  material  (Appendix  D) , 0 
dimensionless  stream  function,  0 
angular  velocity  of  viscometer  bob,  radians/sec 


In  accordance  with  latter  fron  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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